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Abstract 

On a compact Riemannian manifold {Vm,g), we consider the second order pos- 
itive operator = eAg + (6, V) + c, where — is the Laplace-Beltrami operator 
and 6 is a Morse-Smale (MS) field, e a small parameter. We study the measures 
which are the limits of the normalized first eigenfunctions of as e goes to the 
zero. 

In the case of a general MS field b, such a limit measures is the sum of a lin- 
ear combination of Dirac measures located at the singular point of b and a linear 
combination of measures supported by the limit cycles of b. 

When 6 is a MS-gradient vector field, we use a Blow-up analysis to determine 
how the sequence concentrates on the critical point set. We prove that the set of 
critical points that a critical point belongs to the support of a limit measure only 
if the Topological Pressure defined by a variational problem (see |22) ) is achieved 
there. Also if a sequence converges to a measure in such a way that every critical 
points is a limit point of global maxima of the eigenfunction, then we can compute 
the weight of a limit measure. This result provides a link between the limits of the 
first eigenvalues and the associated eigenfunctions . We give an interpretation of 
this result in term of the movement of a Brownian particle driven by a field and 
subjected to a potential well, in the small noise limit. 
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1 Introduction 

Let {Vm,g) denote a compact Riemannian manifold of dimension m>2, with no bound- 
ary. —Ag denotes the associated Laplace-Beltrami operator . In this paper we consider a 
second order elliptic operator depending on the parameter e > 0, 

L, = eAg + e{b) + c, (1) 

where b denotes a C°°vector field on V, c a strictly positive C°° function on Vm and 6{b) 
the Lie derivative operator associated to b: 6{b)u = du{b) for any function u on V . 

is a positive operator to which the Krein-Rutman [25J theorem can be applied. 
Hence the smallest eigenvalue of is simple and strictly positive. The associated 
eigenspace is generated by a strictly positive function m^, normalized in L2(Kn)- The 
behavior of as e goes to zero has been extensively studied [201 !2I1 !221 !S1 !11 • 

In the case when the w-set is a disjoint union of compact invariant hyperbolic set, 
under a mild additional assumption Y. Kifer has proved that the limit of A^ as e goes to 
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zero is the topological pressure (TP), associated to the flow of b and the potential c (for 
deflnitions and details see [22] )• 

Unfortunately, much less is known about the weak limits of the eigenfunctions u^, as 
e goes to zero. By weak limits, we mean the limits of the measures uldvolg in the weak 
topology of measures. Past results include the case when is a bounded domain of M", 
the limit sets of the vector fleld h is reduced to a single point, and are related to the 
zero Dirichlet boundary condition. 

When this singular point is attractive, as e goes to zero, converges to a Dirac 
measure supported by the attractor. When it is repulsive converges to a constant on 
every compact sub-domain of the domain of deflnition (see |21E])- When the attractive set 
consists of a single limit cycle, the flrst eigenvalue converges to the average of the potential 

c of the PDE, along the cycle (see jTHj ), lime^o-^e = ifl-i^^i^^^ll^^ where Xq parametrizes 
the limit cycle of period T, the proof uses a stochastic approach. It is interesting to note 
that this approach and the deterministic methods give complementary results. 
When h = 0,the equation ^ reduces to: 

eAgUs + CUe = XeUe 

The limits of the uldvolg (normalized by: fy uldvolg = 1) has been studied for example in 
([SHEOIEI) and is known as the semi-classical limit. In particular when the potential 
c has a double well in Pi and P2 ( that is an absolute non-degenerate minimum at each 
points), it is well known that the limits of the measures uldvolg as e goes to can 
concentrate on those points only and in the distribution sense, as e goes to 

uldvolg ci5p^ + C25p2, 

where ci + C2 = 1. As far as we know, nobody has addressed the question of computing 
explicitly the coefficients Ci and C2 in general. Of course, under additional assumptions, 
for example, when the are invariant by a group of isometrics, Ci = C2 = |. It is not 
clear at present time whether or not those coefficients are unique. We will see here that 
those coefficients depend only on the Hessian of the potential at the points Pi, P2 and on 
the limit of the ratio (called the modulating ratio) of the local maximum to the global 
maximum of the m^, but they do not depend on anything else. 

Sometimes, as explained in one of the coefficient is zero and the limits of the 
measures uldvolg concentrate on the remaining point. In that case, we use the terminology 
of and say that the degeneracy is removed. The degeneracy can be removed by looking 
at the expansion of in terms of e in the following way. A^ depends on the values of 
the potential and of the derivatives of the metric tensor at the minimum points. The 
Taylor expansions of A^ have to be the same to each order at Pi and P2- If this is not 
the case then only one point will be charged. In this paper we compute the necessary 
conditions up to order 4, in order that the degeneracy can be removed in terms of the 
geometry. These conditions narrow down the set of points which can be charged ( i.e. the 
set of points where the coefficient Ck is strictly positive). 
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We conjecture that the coefficients are unique and when the degeneracy cannot be 
removed they are charged by a global maximum sequence (the modulating ratio is equal 
to 1). 

When 6 is a gradient b = V0, the limits of the measures u^dvolg as e goes to can 
concentrate on the critical points of h only. The concentration results are proved for the 
normalized measures 



e '^^u'^dvolg 
L e~^uldvol 



(2) 



a 



and we shall see how the potential c interacts with the field h to select the points that 
will be charged. This set is the set of points where the Topological Pressure attains its 
minimum (see 

When 6 is a not necessarily the gradient of function, very few results are known 
about the behavior of the measure u^dvolg as e goes to zero, because it not clear by 
which function if any, the function should be replaced in expression |21 Recall that the 
operator with general drifts have not been considered before, especially in the context 
of Quantum Mechanics, because it has an interpretation only in the context of diffusion. 
Actually, on the striking results in this paper is that can be replaced by any global 
Lyapunov function L associated to the field b. Using a weight which is a Gaussian in the 
Lyapunov function, is a crucial input into the problem because it enables us to "filter" the 
limits in order to get results about the concentration. Using a Lyapunov type function, as 
presented in section 01 is by far more general than using a solution of the Hamilton- Jacobi 
equation, 

|VLp + (6,VL) = 0, (3) 

as it is currently used in the formal expansion of the WKB theory. Another striking result 
is that the possible limit measures are supported by specific limit cycles of the field, which 
are not necessarily the attractors.On a compact Riemannian manifold, the existence of a 
solution to the Hamilton- Jacobi equation is not always guarrante to be smooth, while in 
appendix 2, we prove the existence of a smooth Lyapunov function. In this context, we 
prove that the measure 

e~~utdvol„ 

(4) 



Jv^^e '^^uldvolg 

can concentrate on the limit cycle of a Morse-Smale dynamical system, without any addi- 
tional assumptions. Finally, we give explicit results about the decay of the eigenfunction 
sequence near the concentration set and analyze the influence of the Riemannian geometry. 

Our main results in this paper are: 

• In theorem [7| of section 01 we prove that the limits of the normalized eigenfunctions 
Me as e tends to zero, are measures concentrated on the limit sets of a Morse-Smale 
field h. The possible limit measures are supported by specific limit cycles of the 
field, which are not necessarily the attractors. 
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• In theorem |21 we prove that the blown up function is the standard solution of the 
Harmonic Oscillator. This lead us to theorem |21 to give a precise expression of the 
coefficient of the limit measures and to characterize the set of minimum of c which 
are in the support of the limit measures. As a byproduct of the analysis, we obtain 
some estimates of the velocity at which the sequence of local maximum converges 
to a critical point of the potential. 

• In theorenJU we give a geometric interpretation of the coefficients of the expansion 
of in terms of the power of ^/e. These coefficients are invariant of the couple 
(Riemannian metric, potential). The value of the second term of the expansion of 
Xe in powers of ^/e is computed, using minimax procedures, but this value has to 
be compared to the one obtained in |I2|, based on the WKB formula. The result 
obtained here are based on the variational approach and seems to lead to results 
that have to be compare with the results obtained by the formal WKB expansion 
of the first eigenf unction. 

• In theorem we study the case where the field b is a gradient of a Morse-Smale 
function. We prove that the concentration of the first eigenfunction occurs at the 
critical points of h where the Topological Pressure is attained. In addition, the 
weights of the limit measures is given under specific assumptions. 

Remarks. 

Some of the results presented here complete and extend also some previous work 
in analysis [T71 HHJ 12 El El- In addition, at the time the results of this paper were 
announced we were not aware of any reports about the fine selection by the limit of 
the eigenfunction sequence of a subset of critical points of the potential. The detail of 
this selection process is presented in the subsection 12.6.21 and 12.71 and should clarify how 
the metric and the potential are involved. This results extend in particular the work of 

1.1 Notations 

(x^,...,a;™'):?7 — > M, is a coordinate patch on V and / a function on V: 
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1=1 



dx^ 
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1.2 Normal coordinates 



For each PGCminjWe choose a normal coordinate system (xi,...Xjn):Up — > R, centered at 
P, defined on a domain Up such that: 

1. xiX...xxto(C/p) contains the closed ball Bp(S) centered at P and having radius S > 0. 

2. for all i,j, 1< i,j< m, = K{P)Sij. 

3. UpnUg = for all P,Qe C^^,P^Q. 

In the following we will identifie Up with the open neighborhood of O in W^,xi x ... xx^(C/p) 
For r such that Bp(r)cUp,Bp(r) will denote both the geodesic ball centered at P and 
of radius r in V or its image by the mapping xiX...xxm- For the sake of streamlining 
the notations, we shall commit the abuse, when working with the coordinate system 
{xi,...Xm)-Up — > M, of denoting by -Bp(r) the ball of center O and radius r whatever the 
value of r(> of course). 

Now a few words about the blow-up procedure. On a manifold V let P be a point 
and a chart {U,xi,...Xm) of V centered at P: Xi{P)=0, 1 <i< m. The blow-up of power 
t >0 associated to P and the chart {U,xi,...Xm) is the diffeomorphism Bit '■ U R"*, 
Bh{Q)—{\^i{Q) , \xm{Q))- All functions, tensors, differential operators can be trans- 
ported to the open subset Blt{U). Suitably normalized by a power of t, they will have 
limits when t goes to which will be defined on M™. These limits contain a trove of infor- 
mation about the behaviour of the original objects in the neighborhood (more precisely in 
the infinitesimal neighborhood) of P. To simphfy the notations we shall write: |(5 instead 
of Blt{Q), J A instead of Blt{A) if ^4 is a subset of U and so on. 

Recall that a field b is Morse-Smale MS if : (i)the recurrent set of b consists of a 
finite number of hyperbolic points and periodic orbits (ii) each pair of stable or unstable 
manifolds of these points or orbits intersect transversally. If moreover 6 is a gradient of a 
function with respect to the metric g, b will be called a MS gradient field. 

1.3 The self-adjoint case 

In the self adjoint case the vector field b is zero. This assumption simplifies the problem 
because it can be handled by variational methods. Theorem 1 below seems well-known 
but we could not find a proof for Riemannian manifolds, using deterministic techniques 
in the literature. Hence as a starting point, we provide a simple one here. 
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Notations. 



For each PGCmm,we choose a normal coordinate system {xi,...Xm)'-Up — > M , centered 
at P, defined on a domain Up such that: 

(Up) contains the closed ball Bp (5) centered at P and having radius 6 > 0. 

2. for all i,j, 1< ij< m, ^(P) = A,(P)5,,. 

3. UpnUQ = for all P,Q e C^in,P ^ Q- 

In the following we will identifie Up with the open neighborhood of O in M™, xiX...xxm{Up). 
For r such that Bp(t)c Up, Bp{i) will denote both the geodesic ball centered at P and 
of radius r in V or its image by the mapping xiX...xxm- For the sake of streamlining 
the notations, we shall commit the abuse while working with the coordinate system (xi, 
...Xm)'-Up — > M, of denoting by Bp{i) the ball of center O and radius r whatever the 
value of r(> of course). 

Theorem 1 Consider the first eigenvalue problem for the operator Ag + c where c is a 
function with a finite set of minimum points, Cmin, which are not degenerate (in the sense 
of Morse). Assume that the first eigenvalue of the operator is positive. Xs has the 
following variational expression: 



Iv 

X. = inf 



e| I Vm| + cu^ 



dvolg 



u(^m{v)-{o} fyu'^dvolg 
Then, when e converges to zero, X^ converges to the minimum of the function c and the 

u'^dvolg 
fyU^ dvolg 



set of weak limits, when e goes to zero, of the family of measures "^^g defined by the 



positive solutions of the PDE, 

e/S.gUe + cUe = X^Ue on V (5) 

is contained in the simplex 

M = {u = J2hpSp\ Peani„)||5^7P = l,7P>0} (6) 

of all probability measures with support in the finite set Cmin where 6p denotes the Dirac 
measure at the point P. 



Remarks, is uniquely defined up to a multiplicative constant by the Krein-Rutman 
theorem (see P5]). 

In the following proofs, over and over, we will chose appropriate sub-sequences of {ue, 
e > 0}, {Xi;\e > 0} and so on, without saying so explicitly: In order to keep the notations 
simple we will write u^, X^, ... instead of a sequence {u^Jk = 1..), {Xf:Jk = 1..).. 
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Proof: without restricting the generahty we can assume that c > 0. To start note 
that Ag > mine > 0. Using a constant as test function in the functional that define the 

first eigenvalue, we get: 

fy c dvolg 
K < , .... < sup c. 

VOlg{V) V 

Ag is a decreasing function of e because the functional decreases with e. being bounded, 
there exists a sub-sequence which converges weakly and 

V0 e H\V) , / [e(A0)ii, + c(f)Ue] = A, / Ue(f) 
Jv Jv 

To obtain an upper estimate for A^, consider the following radial function defined on 
for m > 2: 

^(m-2)/2 ^(m-2)/2 
= (^2 + ^2)(m-2)/2 - (52 + ^2)(m-2)/2' ^^(^) 



= , on y - Sp((5) (7) 
Then belongs to H^iV). We use this function in the energy function I{u) = '^^^^ '^'^l^u'^a-^i 



jy U-WVUig 

The following standard computation gives an estimate of 



e / llV^^llJcivo/g = e{u^_i [ r'^ ^dr ^ + o(/x)) = ec(m) + o(/i)e 

where the constant c{m) depends only on the dimension n and uJm-i is the volume of 
the unit (n-l)-sphere of R"*, r = d{P, Q) is the geodesic distance and 5 is less than the 
injectivity radius. We have to evaluate the other quantities in the functional. Using the 
change of variable x — y/i, we get: 

/ c4>1^dvolg = / a(j)ldvolg = jJ^V {P)uJm-iJ {5 / jj) + o(/x). 



where 



We have J{5/^l) = for n > 5 and for n=4, 7(5///) = 0(ln5//i) 

then after some computations: 

/ (jpjj^dvolg = f/u;m-iJ{S/ijL) + o{n). 
Jv 



Taking /i = e^^^, we get the expansion 

/(0,) = c(P) + o(l), 

where P is a minimal point of c. 

In the 2-dimensional case, we can consider the test function 

= 0, on V ~ Bp{6) 

where p satisfies < p < 1/2. The same computations as before, taking /i^ = e, gives the 
result. Finally 

mine < A, < c(P) + o(l) 
lim Ae = mine 



Using the energy equation, we have A^ > e |[|Vue||^ + rninc], which forces e jy ||VMe||^ 
dvolg to tend to zero as e goes to zero: 

lime / \\'VuJfdvoL = (8) 
e->o Jy 9 y 

Also for any (f) gC^(V) multiplying equation ^ by the function (j)u^ and integrating by 
part gives : 

[ e||VM,||j0 + CM^0 + e(A0)^ dvolg = X, [ (pu'^dvolg (9) 

Jv I 2 J Jy 

where Jyu"^ = 1- Because the first term and the last term on the left hand-side are 
converging to zero (see|HI), we obtain that 



lim / 0u^(c — \^)dvolg = 0. 
Jv 

< (c — m.mc)(j)u'tdvola < 
Jv V Jv 



But: 

(c — K)(l)nidvolg = 



lim I (c — minc)0M^ = (10) 



e->0 ly V 



Moreover if u denotes a weak limit of as e goes to zero, relation pOj) implies that 
\f(j)eC^{V): 

(c — mmc)(l)u'^ dvolg = (11) 



V 
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Hence if the set {P G V\ c{P) = mine} is of measure 0, we conclude that u = 0. Hence 

all weak limits of are zero and the sequence concentrates, as we shall see shortly. 
Note that relation (jlUj) implies that for any ip G C{V), which is zero in a neighborhood 

of Crviivi . 



£->0 



lim / tjjuldvoL = 



(12) 



If A is a measurable subset of V such that A fl Cmin = 0, A denoting the closure of A, 
applying (fT2|l to a positive continuous function ip with support disjoint from Cmin and 
> 1 on A we get: 



lim / uAvola = 

e->0 ' ^ ^ 



In fact, we can prove that any sequence {en\n G N} converging to 0, contains a subse- 
quence {£nk\k G N} such that the corresponding u^^^ converges to a convex sum of Dirac 
distributions located at the minimum points of c.We assume the u^^ normalized so that 
/y ~ ^- Consider the following decomposition 



buldvolg= V / {{(f)-(f){P)) + ^{P))u^^dvolg+ (pu^dvolg (13) 



Relation ()13|1 implies that: 



[ (t)uldvolg-(t){P) V / uldvolg = V / {(l)-(t){P))uldvolg+ [ 
Jv Pectin -^^p^^) ppr7.„,„ JBp{5) Jv 



(pu^dvoL 



111 



'V-Upsc^.^^Bp(S) 

(14) 

By the continuity of (f), given an 77 > 0, one can find a S{ri) > such that \(f){x) — (f){P)\ < t] 
if xG Bp{6{ri)), for all i. Hence with N=card Cmin: 



Pectin "^^^W")) 



< Nr] 



Relation p2|) implies that: 



lim 

e->0 



bu^dvoL = 



V-Upec,,,i,,Bp(5(rj)) 



After choosing a subsequence of {en\n G N}, still called {en\n G N}, if necessary, we can 
assume that all the limits lim L,<i>u^ dvoL , lim f„ .-.u^ dvoL, P G Cmin, exist. 

Relation p4|) implies that after choosing a subsequence of: 



lim 

n— >oo 



/ (pu^ dvolg — (f){P) lim / dvol 



< Nr]. 



(15) 
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Now note that lim /^^^^^ does not depend on 5. Let ^i, ^2, < 82- Then 



n— >oo 



u1 dvoL = u1 dvoL + / 

Bp{52) JBp{Si) JBp{S2)-Bp{Si) 



Relation ()12p imphes that: 



?i— >oo 



hm / u^^dvolg = 



Bp{52)~Bp{Si) 



Hence if one of the hmits hm f„ ^ u"^ dvoL, hm f„ . dvoL exists, so does the 
other and is equal to it. Set: 



'BpiSiv)) 

Relation p5|) implies that for any 77 > 0: 



lim / ul dvolg = 7p^ 



< Nr] 



lim j (f)u\^dvolg = (j){P)'~fp (16a) 

(j) G C{y) being arbitrary, ()16ap shows that the sequence of measures u^^volg converges 
weakly (in the measure sense) to the measure 'Yl,p^c ■ Ip^p- Finally note that if we apply 
(jl6a|) to the constant function 1, we get: 

j2 ip = i 

111 

Now we prove that supu^ diverges to infinity. Because Jy tends to zero, if supw^ were 

V V 

bounded, then using the following inequality, we get a contradiction with 

1 = / < suptie / u. 



V Jv 



as the right-hand side would converge to zero. Let be a maximum point of u^. Because 
the manifold is compact, it is possible to find a subsequence of which converges to a 
point P. centered at the point P. fe(x) = "'=(v^'^+-^^) and denotes the rescaled metric. 



V 



then the At a maximum point P^, using the maximum principle, c{P^) < A^. Since 
converges to the minimum of v, at the limit, c(P) < mine. This proves that P is a 
minimum point. Using the fact that P is a nondegenerate minimum point, d{P^, P) < 

We will use the following definition. 
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Definition 1 The coefficient 7p is given by, for all S small enough, 



n— >oo 



7p = lim / u^^dvolg 

lBp{5) 



is called the concentration coefficient or the weight of the limit measure at point P G Cmin- 

The coefficient depends on the subsequence, but not on 5. This coefficient characterizes 
the concentration measure at point P. 

1.4 The first eigenvalue problem for the gradient case 

We consider the hmits of the ffist eigenfunctions as e goes to zero when the vector field 
h is the gradient of a Morse function. We estabhsh a result similar to the one obtained 
in the last paragraph : when e converges to zero the limits of first eigenfunctions in the 
weak topology of measures concentrate at the critical points of the field h. 

Consider a Morse function and the vector field, h = V0 and a function c chosen such 
that the eigenvalue Ae of the operator eAg+ < 6, V. > +c is positive on the manifold. To 
study the solutions of the PDE 

eAgU^+ < b, > +cu^ = X^u^,oia V (17) 

we use the transformation b = V0 = — 2eVln'?/'e (it is defined up to a constant) and 
consider the new variable = u^ipe- Equation (fT7j) is transformed into the following PDE 
where the first order term disappeared. 

e'^AgV^ + CgWe = eAe^e, on V 

where = ce + ^ + 

Using the theorem of the preliminary section, we obtain the following results : 

Proposition 1 Suppose that the following condition is satisfied: at the critical points P 
of the function (p, c{P) + A0(P)/2 > 0. Let be a minimizer of the following variational 
problem 

e''Iv[W\l + c^^']dvolg 
eX, = mf — f- — 

v€m{v^)-{o} jyV^dvolg 



then lim^^o ^-^e = niin| | V^l |^ = 0. The weak limits of the normalized measures j ^-^/e^2^^gi 
have their support in the set of critical points of (p. supt;^ tends to +oo as e goes to zero. 

V 

Proof. The proof is very similar to the proof of Theorem Considering the one 
parameter family of eigenfunctions, we obtain that: 

lim^^^ , ' =min V0 ' = 

6^0 J^e-'f'/'uldvolg V ^ 
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and for all function ip G C{V) we have 



where the concentration coefficient 7p. is now defined by 

!bp.(5) e-'^uldvolg 

7p^ = lim ^ , 

''^^ Jy e~~u'^dvolg 

(the limit is independent of 5). The measure ^-4,%^^idvoi converges weakly to YlT=i '^f^Pi 
where Pi are the critical points of the function or the zeros of the vector field b = V0. 
The proof follows exactly the same steps of the previous theorem. | 

Remark. 

The gradient case teaches two things: one is that the concentration occurs on some specific 
sets, related to the vector field and not to c and second that the role of c is to select the 
subset of concentration. 



1.5 The radial case : an example 

In this paragraph we give an example where the recurrent sets of the vector field consists 
of a limit cycle and the sequence of eigenf unction concentrates along this limit cycle. In 
fact, we obtain in presence of radial symmetry a uniform distribution for the limit. 

Consider an annulus A of (A = {x e W\ 1/2 < ||x||Mn < 3/2}), and the radial 
function u^, solution of the partial differential equation 

eAgU^+ < b, > H-CMe = K^e, on A (18) 

Me = 0, on the boundary dA, 

where the field b is given by : 

br = {l-r) 

be = 1 (19) 

and the function a is radial and positive. The field b has an attractive limit cycle at r = 1. 
The problem reduces to: 

u 

e{—drrUe ^—^)+br-drUe + ttU^ = XeUe, OU A (20) 

r 

Me = 0, on dA 

br is the gradient of a function of r. Hence the results of the previous paragraph can be 
applied here. The presence of a boundary does not invalidate these results because the 
limit cycle is an attractor. As e tends to zero, tends to a hmit entirely supported by 
the limit cycle (see also Friedman [2j). 
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2 Blow up analysis with no vector fields 



In the next sections, the hmit measures are analyzed using a blow-up procedure. We shall 
prove that as e goes to 0, the eigenfunctions blow up in the neighborhood of some points 
that are determined by the potential c and the vector field b. The speed with which these 
eigenfunctions blow up can also be determined when b=0 or when b is a gradient field, 
using the Lyapunov functions associated to the field. 

The results differ substantially in the two cases. It appears that the correct scaling is 
not the same in the case when there is only a potential c and the case where there are a 
potential c and a vector field h. 

The general case, where the field can have recurrent sets of integer dimension n > 1 
will be considered elsewhere. The concentration phenomenon is much more complicated, 
depending on the set and on the chose of Lyapunov function. More important, it cannot 
be studied by variational techniques, see [Hj. 

In this section we determine exactly all the possible limits of the eigenfunctions as 
e tends to 0, when there is no field. The main result says that the limit measure is 
concentrated on a subset of the minimum point of the potential c. This limit set is useful 
in the study in the small noise limit, the movement of a random particle moving on a 
Riemannian manifold in the presence of a killing potential c ^3]. 

Moreover, we can explain the assumption 4, p. 93 made by B. Simon ([30]) to study the 
double-well potential problem when e is small. The blow-up method provides a method 
for the explicit computation of the concentration near a bottom well. This generalizes 
also the results obtained in part 9 of 5J about the concentration of the eigenfunctions in 
the case of W^. 

Let us recall the eigenfunction problem, 

eAgU^ + CUe = XeUe (21) 

/ u'^dvolg = 1, 

c is a Morse function and Cmin, denotes the subset of minimal points. Recall the quotient 

g ^ Iv [g||V^llg + CM'] dvolg 

fy u^dvolg 

We shall now state and prove the main theorems of section El We introduce some concepts 
which will be used in the proofs of these theorems. 

Now a few words about the blow-up procedure. For each P G Cmin; we choose a 
normal coordinate system (xi, ...Xm)'-Up — > M , centered at P, defined on a domain Up 
such that: 

1. xi^...^Xra{Up) contains the closed ball Bp{5) centered at P and having radius 
6>0. 

2. for all i, j, 1< i, j< m, ^(P) = A,,(P)<5,,. 
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3. C/p n C/q = for all P, Q e C^^, P^Q. 

In the following we will identifie U p with the open neighborhood of in M™, xi x...xxm{U p). 
For r such that Sp(r)c Up, -Bp(r) will denote both the geodesic ball centered at P and 
of radius r in y or its image by the mapping xiX...xXm- For the sake of streamlining 
the notations, we shall commit the abuse while working with the coordinate system {xi, 
...Xm)'-Up — > M, of denoting by Bp{y) the ball of center O and radius r whatever the 
value of r(> of course). On a manifold V let P be a point and a chart {U, xi, ...Xm) of V 
centered at P: Xi{P)—0, 1 < i < m. The blow up of power t >0 associated to P and the 
chart {U, Xi, ...Xm) is the diffeomorphism Bit '■ U — >• M"*, Blt{Q)={jXi{Q), ...,jXjn{Q))- 
All functions, tensors, differential operators can the be transported to the open subset 
Blt{U). Suitably normalized by a power of t, they will have limits when t goes to which 
will be defined on M™. These limits contain a trove of information about the behaviour 
of the original objects in the neighborhood (more precisely in the infinitesimal neighbor- 
hood) of P. To simplify the notations we shall write:|(5 instead of Blt{Q), jA instead of 
Blt{A) if is a subset of U and so on. 



In the following all the blow-ups will be associated to geodesic charts ([/, Xi, ...Xm) 
with pole at P. On the magnified set -^xiX...xxm{Up), we can define the function wp,£ 



blow up of the function^: 



where Up 



-max Ug. 

V 



wp,e{y) 



Up 



2.1 Main theorem 



Definition 2 We define A as 



A = inf 



n=l 



Theorem 2 Selection-Concentration. 



• (i)For any Pe Cmin,fl«?/ sequence of Vs's, with e tending to 0, contains a subse- 
quence {tt£„}; such that the sequence of blown-up functions wp^^^ at P converges to 
a function wpM'^->'M.+, both in the L?norm and the CP° topology. 

• (a) w satisfies the equation and inequality: 

m 

Aew + \i{P)x1w = \w 

i=l 
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< w < maxty < 1 
where Ae is the negative standard Euclidean Laplacian on R"*. 
(iimj:n=iVUP)>^ - then 

wp = 

(iv) If Xl™=i V ^rt(-P) = A and there exists a sequence S C N, such that each 
Uir„,n eS, has a maximum point Qn with the property that the sequence {Qn\n& S} 
converges to P, then this function wp is: 



wp{x) = exp 

n=l 



and 

A = A 

VYl^=i V ^n{P) = A and no such subsequence S exists, then 

A = A 



ni 



n=l 



where fp is a factor > 0, which depends on the sequence {u^„}- 

• (vi) For any sequence e's, there exists at least one P e Cmin o-nd at least one 
subsequence S for which the case (iv) occurs. 

Remark. Note that in case (iv) the hmit w is independent of the sequence {ue^}. 

In order to prove the main theorem, wc need two propositions. The first gives estimates 
of the first eigenvalue and the second, estimates of the decay of the eigenf unctions. 

2.2 Auxiliary propositions 

Proposition 2 The first eigenvalue A^ satisfies the following inequality 

mine < Ag < mine + Ae^^^ 



V V 



where A = inf 
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Proof. For all ueH^(V), Qe{u) > min c f^u^. Hence mine < A^. To prove the right 

V y 

hand-side inequality, we will use a test function in the variational quotient Q. In the 
neighborhood of a point P e Cmin, consider the function: 

(p^ = e~^''*^</^ - e~^, on Mp{p) 
= 0, on V-Afp{p) 

Mp{p) is the connected component of the set {x|^™ ^ \/\.x} < p} containing 0, where 

p is taken so small that Np{p) is contained in Bp(5). We take as coefficients pi — ' 
where \ — \{P) for simplicity. In the coordinate system at P, 

9iAP) = Sij + 0{\\x\\l^), 



Ric denotes the Ricci tensor. We denote nrA*i A* ^'^'^ the quadratic form Ylil^ii'^^Y 
by q. We recall that: 



/ 



TT 



X 



dx = . — , 

V fJ'i 



2p. 



3/2 ■ 



To evaluate the quotient Qe{u^), we compute the leading terms in e of the integrals p 



I 
I 



\\V(f)e\\ldvolg 

ccpldvolg, 



P(p) 



! 

Ja 



Afp{p) 



{\\V(t>,\\l)dvolg^ 



E 



Np(p) —1 dxi dxj 



^/detgdx, 



where g''^ is the matrix inverse of gij and dx is the Lebesgue volume. 



/ 

Ja 



J^p(p) 



^ ] l^i Xi + ^ ] OiijkXiXjX}; 

i,j,k=l 



i=l 



e '^dx, 



where the functions aijk are defined C°° and bounded on Bp{6). Performing the blow-up 
at P i.e. the change of variable Zi = Xi^/pi, 



f {\\VMl)dvolg= [ 
JNp{p) Jb 



B(0,^) 



Zi Zj Z]^ 



(II 1 1 2 \ dz 
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where // = UlLi l^i- 



e [ {\m\\l)dvol,^^\±^ + 0{ei). 



We will now evaluate the potential term: 



/ ccfy^dvolg — j 
JMp(o) J a 



Afpip) 



iP) + X,x^ + 0{\\x\ II.) I (^1 - ^^^^xV + 0{ 



X 



(P) / <t>ldx+ f 

JUp{p) J A 



Mp{p) 



.k=l 



c{P)Ricu{P) 
6 



i,.i,k=l 



dx 



dx 



where the functions a,ijk are defined C°° and bounded on B{S). With the same change of 
variables, after expanding the square = — 2e~*/^~''/^^ + e~^/^ we get: 



/ 



x1(f)1dx 



m/2 



2^V Ai(P) 



Otu(exp— ^) means that for any 77 e [0,1[, there exists a constant K{r]) independant 



of e sucht that the error is at most equal in absolute value to K(77)exp 
symmetry argument shows that: 



/ XiXj(f)'^dx = if i 7^ j 

JAfpip) 

The potential term in the integral becomes: 

JAp(p) yH' \ 

Also: 



VP 



. An easy 



JAfpip) V/^ 



1=1 



P III 



Ricii{P) fT 
12~V A" 



6 Vr+«(^ 



The quotient can now be evaluated: 



(Er=i^+c(p)+Er=i(A.-c(p) 



^) /T + 0(ef) 
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c(p) i-Er=i^^/i +Er=iv^+o(ei)) 

Qe{4>e) = — 



If we use the fact that c(P) = mine, the quotient can be simphfied as follows 



V 



If we take the minimum over all test functions centered at any critical points of the set 
of minimal points Cm, we obtain the estimate: 

mine < A, < < mine + inf J ^ ^\i{P)t\P G C„,^A+0{e) < mine+ Ae^/2 + 0(e). 



i=l 



Remark: We expect in general for smooth potential that there exists an asymptotic 
expansion: 

k=0 

If it does exist can one find a systematic procedure to compute the coefficients c^? From 
the previous result, we have that 

Co = min e, 

which is the Topological Pressure. We will see in the following results that 



Cl 



1=1 



Proposition |21 provides an estimate of the velocity of convergence of the sequence of 
maximum points of u^ to an element of Cmin- 



Lemma 1 (i) If for a sequence e„ tending to zero, lim„^oo jy'^l^'^'^^h > , t/i 



en 



lim supMe„ = +00. 

n^oo y 



Proof. (i)Suppose that for a subsequence e„, still denoted by e„, lim„^oo supy u^^ < +oo. 
Then for all n, supy u^^ "£ N, a constant. For any 77 > 0, chose an open neighborhood K 
of Cmin, such that volg{K) is smaller than Now 



dvoL = dvoL + dvoL 

V JK JV-K 
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I u^^dvolg = / u^^dvolg + snpUe^ j — — — dvolg 



By Appendix II, ^^^"^ — >0 uniformly on V — K. Hence limn^oo JyU^^dvolg < t]. 

Because rj is arbitrary, lim„^oo Jv'^'i„^'^'^^9 = 0. A contradiction. 

Recall that dg denotes the Riemmanian distance associated to the metric g, 

Lemma 2 For each e, let us denote by Aie the set of all maximum points of v^. There 
exists a constant A depending only on c, such that: 

sup dl{q,C^in) < Ae^^^ 

It follows from this that the set of limit points of the set of maximum points of is 
contained in CLin- 



Proof. Let qG Ais- Recall eAgU^{q) + c{q)u^{q) = X^u^{q). Because the solution u^ is 
positive and AgUe{q) > 0, at the maximum point q of Us, c{q) < Ag. By Proposition 121 
< c(g) — mine < — mine < Ae^^^. Because the critical points of c are non-degenerate, 

it is easy to see that there exists a constant F depending only on e such that for P E V, 
d^(P, C„nn)< r(e(P)-min e). Take v4=rA. | 

Remark. This result proves that any sequence of e's converging to 0, contains a subse- 
quence {ek\k e N} such that there exists a P G Cmin and a vector P*G TpV with the 
property: 

P, = expp(ei/^P* + o(e^/^)). 

The length ||P*|| of the vector P* is the distance between the peak of concentration and 
the set Cmin in the blow up space. It can be considered as a measure of the convergence 
velocity. 

We have so far computed an estimate of the rescaled eigenvalue ^iz£H^IliiL£_ Now we 
will provide an estimate of the eigenfunction in the neighborhood of the points in Cmin- 
Let P be a point in Cmin- Recall that = and = max u^. 

Proposition 3 For all eq g]0, 1[: 

• (i) sup Jbp(s/ ^)'^'^^y^^^y ^ +00- More generally, for any continuous function 

]0,eo] 

f : [0, 1] X R™ — >M, (e, — >f{e,y), having at most polynomial growth at infinity, 

sup/ f{e,y)we{y)'^dy < +00 

]0,eo] J Bp{S/^ 
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(ii) the set of restrictions wl\Bp{5/^), e g]0, 1], of the to the halls Bp{5/-^) 
satisfies the following condition: for any rj > 0, there exists a compact K C M™ and 
a eijj) > such that 

f{^^y)we{yfdy < r], 



Bp{5/^)-K 



for all e g]0, e{7])\. 



2.3 Proof of the Proposition [HI 

To start with, rewrite equation (j25) a^s follows 



v/eA.u, + -^^ = 



Introducing the notations = ^ ""/iT" and /Xe = for simplicity, we have: 

The function e~*^= "^£(2;) is the solution of the parabolic Cauchy problem: 

^ = -yfe^gP - c,p 

p(0, x) = 

Note that the sequence /i^ is bounded. We estimate in the ball Bp{6), using the fact 
that the restriction of the function e~^'*fe(x) to Bp (5) is the solution of the parabolic 
initial- boundary value problem: 

^ = -V^^gP - C,p 
p(0, x) = Ue{x) 
P{t,x)ldBp(5) = U,{x)ldBp{S)e~'''^. 

In the coordinate system at P, 



dt ^ dxidxj ^ dxk 

i,j=l •' k=l 

where 

m 
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For xG Bp{S), the solution is given by the Feynman-Kac formula with boundary term: 



e ^'^uJx) 



(22) 



where ^^(t) is the process starting at x at time and satisfying the Ito equation: 

dX,{t) = y/^B{Xe{t)) + ^a{Xe{t))dW{t) for t < rf (23) 

where W{t) is a standard m dimensional Brownian motion and a:U — s>End(M™') is the 
positive definite square root of the matrix function {2g^^). is the first exit time from 
the ball Bp{6), of the process X^, starting at x. Let Me = max u^. Then equation (j^^ 

implies the inequality: 



(24) 

for X G Bp{5). Let us define: 

We estimate the terms I, II independently. 
Estimate of / 

To estimate the boundary term I in inequality (j2H), we apply the results proved in Ap- 
pendix 2 to the equation()21|) taking 'il)=c-vamc, 6=0, q = min c-A^ and e = e. Then for 

any integer k, any compact subset C disjoint from the set Cmin, there exists a positive 
constant A(k, C) such that for e g]0, 1]: 



maXMe < ^(^; C)^ 

Taking C =V-U{Bp{6)\P G C^n,} : 
Estimate of II 

Let a be a number in ]0, |[. Assume that e<l. We split II as follows: 

II = E, (xiX(t<.,^)e-^^'(^^(^))'^^) + E, ((1 - Xi)X(*<rf)e-^o '^^^^^^^^^'^ 
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where Xi is the characteristic function of the set <| sup||Xe(s)||iRm > Se" 

[0,t] 



II = III + IV 

Because c^ is non negative and a < |: 



III = E,., 



sup \\Xs{t)\\Rm. > fe^ 



We have: 



X,{t A = / " B{Xe{s)) + ^ / a{Xe{s))dW{s) 
Jo Jo 



Now sup < 5. Hence ||v/^/o ' 5(Xe(s))||Km < tMiy^mT, where Mi = 

[0,tATf] 

sup ||-B(x)||Km. Then a well known lemma (see [34j) shows that if 6e" > ||x||iRm + 

xe-B(<5) 

tMl^/me, 



P. 



sup ||Xe(t)||jj„ > fe" 

[0,tAT^] 



< 2mexp 



[Se" — \\x\\Rm — tMiy^me 



2mMH^ 



where M= sup ||(^(^)||m^* It is clear that: 

Bp(S) 



Hence if | |a;| Ir™ + tMl^/e < 5e": 



2mMHy/e 



To estimate the second integral IV = E^ ^(1 — Xi)X{t<Tf )e -^o cj(Xe(s))ds j ^ fjggne the 

process Y^{t) for t < = first exit time of the process X^ from t/p: 



Then: 



Y,{t) = YM + ^ / + / a{X,{s))dW{s 
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for t < .But using the definition of a and the system of coordinates (xi, . 

a{x) = V2{Idm + 
where $:C/ — >End{W") is a C°°matrix function such that : 

m 

= XI ^ijkl{x)XkXi 
k,l=l 

Its value at P is: 

^ijkl{P) — ~-^{Riklj{P) + Rilkj{P)), 

where F^) is a matrix function such that: 

m 

%{e,YM) = J2 ^m{</^Y,{t))Y,,k{m,i{t) 

k,l=l 

Also the components of the field B : 

m 

= ^Bij{x)xj 
j=i 

Hence: 

B{X,{t)) = ^B{e,Y,{t)) 

m 

Bi{^,y) = ^Bij{y^)yj 

Finally: 

Y,{t)^Y,{0) + V2W{t) + Z,{t) 

Z,{t) = / B{e, Y,{s))ds + V2i / $(£, Y,{s))dW{s) 
Jo Jo 

To estimate IV, we split it into V and VI, choosing a /3 e[0, a[: 

Y^E, ( (l-xi)X(t<r,^)e--^o-^(^e(^))<i^;sup||Z,(s)||M^ ) 
V [0.*] / 



[o,t] 
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Estimate of V 

We shall use a variant of the inequality (37.9) stated in (|2H1 p.78): 

P{sup ||M,|| > y- [M]t < K} < 2mexp {-^^\ , (25) 
[o,i] V 2mify 

where M is a continuous local martingale, at t=0. We apply this formula to Mt = 
Ze{t) - ^ jQB{e,Y,{s))ds = ll%{e,Y^{s))dW {s) stopped at r^^, the exit time of 

the process from Bp{S). 
Note that : 

(l-Xl)||v^$(£,n(s))|Und(R'") <C25V° 

and that: 

(1 -xi)||v^5(£,i;(s))||m™ < C3fe-+^ 

for some constants C2 depending only on the values of a on -Bp (5), C3 depending only 
on the values of the vector field B on Bp{6). Thus we have: 

{sup ||M,|| > = 0} C {sup ||M,|| > [M]t < t(C25V")2} 

[o,t] [0,t] 

Recall that: 

[M]t = 2^ / tr[$ie, F.(s))$(£, Yeis)y]ds, 
Jo 

[M]t<2 / ||^$(£,n(s))|| 
Jo 



Hence: 



y<PAZe{t)>e^;xi = o} 

V < P{sup ||M,|| >6^ - tC3fe"+^; [M]t < t(C25^£^°)^} 

m 

V < 2mexp-^ ^t,, , ■ 

Estimate of VI 

By Taylor formula there exist a constant Ci > depending only on the function c such 
that for xG Up: 



|c(x) - C(P) - ^\nXl\ < Ci||x||]Rm A^X- 
n=l 1 

Then, recalling that c(P)=minc and Cg= ^~™^^"' ^ for t<T^: 

m m 

(1 - ^Ci||n(t)||M™) ^ A,y;,(t)^ < c,(x,(t)) < (1 + ^Ci||n(t)|k™) 5^ A,n,(t)^ 
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For the sake of simplicity let us denote the positive definite quadratic form "^i-^f 
q{X). Given any a > 0, any a, 6 G R*" we have: 

(1 - a)q{a) + (1 - -)q{b) < q{a + 6) < (1 + a)q{a) + (1 + -)q{h). (26) 
a a 

Now: 

q{Y,(t)) = qiXM + V2W{t) + Z,{t)). 
Using the inequality ^ taking 0=^^(0) + \/2W{t), h=Ze{t), a = I : 

qiUt)) > IqiYeiO) + V2Wit)) - (27) 

Hence setting y=-|= = F£(0),for yG i?(fc""3): 

VI = i5j(l-Xi)x(*<.,^)e-^o'=^(^^^(^»'^^ sup P,(.)||k..<5^ 
Vl<Eyl{l-xi)Xit<r,^)exp -{1-64^0 [ q{Ye{s))ds; sup \\Z,{s)\kn. < eA 

\ JO [OMr^] J 

Vl>Ey({l-Xi)X(t<r^)exp-{l + 6^Ci) [ qiY,{s))ds; sup \\Z,{s)\Wm < 

\ JO [0,iATf] J 

Using the inequality (j2Zl): 

VI < Ey \{1 - Xi)Xit<rr) exp(l - S^^C) [ q{Z,{s)ds- 

Jo 

{l-S^Ol- [ q{y + V2W{t))ds; sup \\Z,{s)\Urn < 
VI < Ci{6,e,t)Eo (^exp-(l - 5^Ci)^^ g(y + y2iy(s))cis^ 

where: 

Ci{S,e,t) = exp ((1 — 5\^Cl)t^/e max{q{v)\ ||f ||ir™ < e^}) 
C4{6,e,t) =exp (^1 - 5^Ci)£5+2/3niax{g(t;)| < 1}) . 

Note that C4{6,e,t) — >1 as e — >0 (keeping t fixed). Let us estimate the expectation 
Eo (exp -(1 - 6^Ci)^ /o q{y + W{t))ds^ . For simphcity set a{e) = (1 - S^Ci). Then: 

Eo (^exp £ q{y + V2W{t))ds^ = U^.E^ (^exp ^ A,(i/, + v^iy,(s))^ds^ 

where y={yi, ...yj and W(t) = (Wi(t), W^(t)). 

To find the value of Eq (exp — ^ /J Xi{yi + V^Wi{s)ydsj we use the following lemma: 



27 



Lemma 3 Consider w{s) the Brownian motion in R, a positive real number A > and 
the function z^: 



where is the expectation for a process starting at a point x. Then is a hounded 
solution of the parabolic equation, 



dz 1 d^z 



dt 2 
z{Q,x) = 1. 



- \{x^ + iJix)z, for X eM.,t>0 



For all X e R,^> 0,the solution is given by 



A/cosn(tv2A) ^ 



and the value of the expectation at the origin is, 



1 



cxp 



cosh(tV2A) 



~ [2) ^^T'' 



V2 



(28) 



Assuming the lemma we pursue the evaluation of : 



Eo (^exp £ K{y, + V2Wi{s)fds^ 



= e ^ Ey,\ exp 



Eq exp 



a{e) 



K{yi + V2Wi{s)yds 



a(s)t\,,,f 



COS 



2V2 



and 

Eq ^exp 
Finally: 

VI < 



ale] 



K{yi + V2Wi{s)fds 



exp 



a{£)Xi tanh(t^ 2a{e)Xi) 2 



2V2 



Vi 



A/ COS h(V2a(£)A0 



C4(5,£,t) 



n-,Jcosli(V2a(£)Ai) 



\J a{e)\i tanh(t^2a(£:)Aj) 2 



2\/2 
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End of the proof of Proposition [31 

Now we can wrap up the proof of Proposition El For any a g]0, |[, any [3 G]0,a[ , any 
e g]0,£:o] wtiere SqKI depends only on a, [3,5 and any integer k, we get tlie estimate if 



Ue n™iVcosh(ty2^) 



exp - ^ 



^/a{e)Xi tanh(tA/2a(e)A 



i=l 



2 



+ 2m 

Define the function We:^p{5) 



Then for yGBp(5/£3-"): 

^.(,)e--* < A{k, 5)e' + ^ ^^f^'^'^ exp (- V^^^^^^^^v/^^^ 



+ 2m 



exp 



n^iVcosh(tV2aAi 

^3^ 



2^2 



+ exp 



2mtM'^y/e 



We want to estimate sup{we(?/)| yG -Bp((5/£:4)}. 



supK(y)|y G fip(5/£3)} < supK(i/)| y G fip(5/^^) " fip(5/2£3'")}, 
supK(y)|y G 5p(5/e^)} < sup{w;e(y)|y G 5p(5/2£^-")}) 



sup{t/;,(l/)| y G Sp((5/£3) - Sp((5/2£^")} < sup{^| x G 5p(5) - i?p(fe72)} 

Me 

Using Appendix 2 as in the evaluation of the boundary integral I, taking e = b=0, 

lb = c — mine, c^ = min c-Ae we get for all n > 1 

^ V V ~ 



sup{4^^| X G Bp{6) - fip(fe°/2)} < C{n 



Me (mm 

where min ip is the minimum of ip on Bp{6) — Bp{6e°' /2) and C(n) is a constant depending 
on the data g, c and on n but not on e. Now there exists a constant Q>0, such that for 
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all X gV, c(x) > mill c+Qdc,(x, Cmin)^. Hence min ijj > Since a < |, we see that 

for any kG N, there exists a constant ■y{k) > 0, depending on 6, such that: 



sup{^^| X G Bp{6) - 5p(fe"/2)} < 7(fc)£* 



Hence: 



sup{w;.(l/)| y G i?p(5/^3) _ i?p(5/2£i-")} < ^{k)e'' 
Finally for xG Bp{6): 



lit, ^cosh(tv/2^ 

+ (A(A:,5)+7(/e))e'= + 2m 
Equivalently, for all yG Bp{S/^): 



I exp I y ^ ■;= I -|- 



i=l 



nl^i V cosh(tA/2a(e)Ai 



/ _ ^ a/ a(g)Aj tanh(tA/2a(£)Ai) j _^ 



i=l 



+ (A(A;,5)+7(A;))e' + 2m 



exp 



2mtC|5%4a 



+ exp 



2V2 



(29) 



2mtM2v^ 

Given a function f as in the statement of Theorem |2l for any integer such that sup 



\fie,y)\ 



< +00, for any integer k > ^ + N inequalitvip^implies that 



sup/ \f{e,y)\w{yfdy < +00 

]0,eo] J Bp{S/^) 



This is the statement (i) of Proposition IHl 



(ii) Fixing t and choosing an integer k>^ + A^, we have: 



\imvol{Bp{6/4^) 

£->0 



1 + 



6 



N 



A{k, 6)e'' + '^{k)e^ + 2m exp 



+2m exp 



[^-tMiy/j)]- 

2mtM^^/e 



Hence there exists an ei{ri) such that for all e g]0, £1(77)]: 
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1. (e/^-tC3fc")2>ip, 



2. [^-tM,^)r> , 

3. a (e) = (1 — 6^Ci) > ^ As a consequence there exists an e{ri, t,k),0 < e{ri, t, k) < 



2' 

Si{r]), such that for e E]0,e{ri,t,k)] 



4. vol{Bp{6/^) \^{A{k,6) +-f{k))e'' + 2m 
1 

Then we can choose a ball B'"(P, R) such that for all e G]0,e(?7, t, A;)] : 



/ 



2 



C4(5,e,t) ^ v/a(£)A.^tanh(tv/2a(£)Aj) ^ \ j 2*^.^ / 



i."^(o,K) y^™^ Jcosh(tv/2^) 



exp 



These two inequalities prove the statement (ii) of Proposition El I 

2.3.1 Proof of Lemma El 

The fact that the function satisfies equation (j28|l is a consequence of Feynman-Kac 
formula (|22|). Making the change of coordinate x — > a; — /i/2, the problem becomes: 

dz 1 d'^z 2 ^ \ 

z{0,x) = 1 

To solve this equation, by the uniqueness of solutions for the Cauchy problem for the 
parabolic equation (|21|)5 it is enough to find a solution of the type e"*^^*-*^' +'^(*). A simple 
computation leads to the two coupled equations: 

+ 202 = A 
^ + (j) = — 

which can be solved easily using the initial condition at time zero, 

0(t) = y^tanh(ty2A) 
^(t) = logcosh(tV2A) + ^ 

so we obtain the expression 

.Vx,t) = ^ ' exp f_ v^^^^Mtv/2A) ra^ , ^ 



cosh(tV2A) 



V2 V 27 4 
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2.4 Proof of the main theorem 



We use the notations of the preceding sections. Because normal 
coordinate system centered at P G Cmin, the eigenfunction Ws satisfies the following 
renormalized equation, derived from equation (PT|) . in the blown up coordinates, y = i^ex: 

c{y^x) — mine — mine ^ ^ 

Ag^w, H ^ We = w„ in — = —xi x ... x x„(f/) (30) 

where is the metric rescaled by e, converging to the Euclidean metric uniformly on 
every compact set of M"* in the C°° topology: 



i,j=l ^ ^ i,j=l ^ 

Remark that all coefficients of the partial differential equation (j3(Jj) are bounded. In 
particular this term is bounded with e. By Lemma El the quantity is bounded 

and nonnegative. 

By the classical theory of elliptic partial differential equations, since is bounded by 
1 for all e, any sequence {we„} with e„ tending to contains a subsequence, still denoted 
by {we„} for simplicity, which converges to a solution of the following elliptic equation 
in the C°° topology: 

m 

Aew + Xnvlw = Xw on (31) 

n=l 

A=lim ^f^^i^ 

n — >oo £ ' 

tn 



where: 



An — A„(P) 



and for simplicity: 

For every compact Kg M"^ : 

/ w{yfdy= lim / WeSvfdy < sn^ / w e{yf dy < +oo 

JK ^—>^Jk ]0,eo] J Bp{5/^) 

Hence wG L^(M'") and I |w| |/^2(]Rm) < lim ||w£„||j;^2(iRm). In fact w= lim We,^ in L^(M™). 

n >oo n >oo 

To see this assume that inf | \ws„ — w\ liamm) > > 0- By (n) of Proposition El and the fact 

n 

that wG L^(M'"), we can find a compact K and an integer Ni such that fj^m_K '^eSuYdy < 
{{f li n> Ni, and J^^_j^ w{yfdy < (|)\ 
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Now there exists an integer N2 such that jj^{we^{y) — w{y))^ dy < (|)^ if n> N2. 
Then for n> max(A''i, A^2), Wwe^ — tf] |^2(]jm) < 3 (|)^ .This contradicts the assumption inf 

\\We^ - U^||l2(R-) > ^ > 0. 

Since w is not zero, equation (j3T|l is an eigenvalue problem. To determine the solutions 
to this problem, consider the unbounded operator L:D — ^ L^(M™), 

m 

Lw = ^ew + \y}w (32) 

i=l 

and 

n=l 

L is a self adjoint operator, the spectrum of which we want to compute. 

Let us introduce the one dimensional unbounded operators L^'-Di — >L^(M) where 
Di={u e L\R)\u G H^{R),ux^ e L\R)}. L^u = -0 + Xny^u, L„ is a self adjoint 
operator, the Hermite operator. It is well known that L^(M™) can be identified to the 
m-fold projective tensor product L^(M)®L^(M)(8)..., (8)L^(R).Then L is the self adjoint 
extension of the operator Yl^=i Ln : Di ® ... ^ Di — >L^(R'"). 

J^L„ {ui{yi) ® ... ®Um{ym)) = ^ui{yi) (g) .Ln{un){y 

vn=l / n=l 

^From this it follows that the spectrum alL) of L is given by the formula: 



a{L) = a{Li) + ... + a{L„ 



where cr(Li) + ... + o-(Lm) is the closure of the set{/ii + ... + Hk G o-{Lk), I < k < m}. 
In the present case: 

a(L„) = {(2fc + l)v^|kGZ+} 

The eigenfunction (up to multiplication by a scalar) corresponding to the eigenvalue 
{2k + l)^/X^ is the function : 

e^^^Hk{y^/K) 
where is the k*'' Hermite polynomial. 

dx'^ 

Hence the lowest eigenvalue of the operator L is An associated eigen- 

function of L is: Y[n=i ^^P (~ j . Up to multiplication by a scalar this is the only 

eigenfunction of L associated to Yl^=i ^'^d it is strictly of one sign. This follows from 
Friedrichs' theorem and the fact that L is a positive self adjoint operator (see |2Z] vol.4 
p.207, Thm. XIII.48). By Proposition [21 : 
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A < A 

Hence we get: 

A = A 

and: 

m 



w[x) = _[ _[ exp ( ^ — 1 = exp - 2^ — ^ — 

n=l ^ ^ ?i=l 

Definition 3 Cminmin is the subset {C\P G Cmin, X]r=i a/ ^n{P) = A}. 

The following theorem amplifies theorem ^ and sums up all our results in the self 
adjoint case. 

Theorem 3 Distribution of limit measures. 

(i) Using the notations of Lemma\^ sup c?^(g, Cminmin) < Ae^^"^ 

(ii) Let S he the weak limit set of the measures u'^dvolg/ Jyuldvolg as e goes to zero. 
Then: 

S = {fie M(\/)|/i = J2 {7P5P| P e C^inmin} , 

(iii) For any sequence {u^,^} such that the measures ul^volg/ Jyul^dvolg converge weakly 
to fi, the concentration coefficient 7p, due to the normalization condition in the space, 
is given 

k=l 

where 

SUPs 

fp = lim 

n->oo snpy^Ue^ 

and 

m 

P&Crnin min ^ — 1 

(iii) If P is a limit point of the sets of maximum points of the Us„ then fp=l. 

(iv) There is always a P such that fp=l 

Proof. We shall use the previous notations. Note that: 

/ uldvolg = 2, / uldvolg + / u^dvolg 



in 



where 

Cmm(^) = Up(zC^.^Bp{6) 
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/^From now on we normalize the requiring that jyu1dvolg=l. Let {ue^|} be any se- 
quence such that the measures v?^^volg/ Jy^^ dvolg converge weakly to /i gM(V). Equation 
imi implies that: 



e->0 



lim / dvoL = 



Hence: 



For any Pe Cn 



^ Pen ■ J Bp{S) 



.^i^n 5Z / uljvolg = 1. (33) 



u'.dvolg = u\e^J^ [ w^^Jvolg^^ , (34) 

Bp{5) JBp{5/^) 



where wp^^ is the restriction of the function to Bp (5) extended by outside Bp{S). 

Any sequence of e's converging to zero contains a subsequence( e„) such that for any 
P in Cmin, lBp{s)'^'^n^'^'^^9 couvcrgcs as n goes to oo and wp,e„ converges to wp{x) = 

fp exp 

(- Er=i ^"^y^"" ) by Theorem 121 Also by the same theorem fp=l for at least one 
P. Hence the sequence {ule':% e n} will converge to a constant K> 0. We claim that: 

lim [ wl^dvolg^^ = [ wp{yfdy = fj 



To see this note that: 



Bp (5/^ JBp{S/^) 

Using Morse's lemma on Up: 

m 

gij{x) = 6ij + ^ Gijki{x)xkXi, 

k,l=l 

m 

Qeijiy) = g-ijiy^fe) = 6ij + ^ G,jkiiy\^)ykyi- 

k,l=l 

Hence: 



det{geij){y) = 1 + y/eG{e,y) 



where: 



\G{e,y)\<C{M){l + \\y\\l^) 
C(M) is a constant depending only on supdGjjfcK^)!!!— < m,xG Bp{S)}. 

/ ^p,eSy)\/^^^((9en)ij)dy = wl^^{y)dy+^/£ wl^^^^{y)G{en,y)dy 

JBp{5/^) ^ JBp{5/^) JBp(5/^) 
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Proposition 121 (i) implies that: 

lim / wj, {y)G{en,y)dy = 



and: 



Hence: 



lim 

n— >oo 



Bp {SI 



Wp.,Sy)dy = fp 11 exp dy 



lim / Wp^ dvoL^ = fp — , 

lim / ul^dvolg = 1 



lBp{5) 

We have then the following alternative for any P E 



mm mm 



f , f o,ifp^ a _ __ ____ 

hm / U^Jvolg = < fl (27r)^ if P c p 

"->-Jbp{5) [ ^^„^^y3;;^,lt^e^mmmm 



Because hm ^p^^^^^^^ /^^^^^ M2di;o/g = 1, 



K= y: 



Pec 

mm min 

Finally LemmsEl implies (iv). | 



Remark 1: As we shall see in the next section, the Topological Pressure gives enough 
information to determine where the concentration process occurs. This information is 
contained in the quantity Cq only. But what we actually found is that a second quantity 
Ci which carries some second order information about the potential, narrows down the 
possible concentration set. 

Remark 2: When c is not a Morse function, it is an open problem to compute all 
the limit and the values of the coefficients. 

Remark 3: We insist that not all minimum points are necessarily charged, except 
when the equations admit some symmetries. It is an open problem to prove that indeed the 
minimum points of the potential c where the Topological Pressure achieves its minimum 
are all charged. 

Remark 4: The main theorem is useful to interpret the dynamics of a particle on a 
compact manifold moving under the influence of small random noise. Indeed , as e goes 
to zero, the density probability of the particle satisfies the Fokker-Plank equation: 

dp 

— = —eAp — cp. 
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where c is a killing term ^H]- P can be expanded in a series of eigenfunctions, 0^ associated 
to the eigenvalues Aj, 

i 

The previous analysis proves that the particle will most likely be found near some specific 
minimal points of c. In general, it is well known by the standard Exit problem (see |29j ) 
that a particle subjected to an attractive field can escape outside the domain of attraction, 
under the influence of small random noise. 

But in the case of a compact manifold, they are several bassins of attraction. The 
particle wanders in and out of these bassins. But, the particle will most likely be found 
in the neighborhood of a point of Cmin, where the topological pressure is achieved. 

Corollary 1 Let W C V an open subset. Assume that dW fl Cmin is empty. If for a se- 
quence e„ tending to zero, lim^^oo [/vi/ '"en^'^'^^ff/ /y '"e„'^'^''^ff] ^ ^' then lim^^oo supw^^ = 

w 

+00. 



Proof. After taking a subsequence if need be, we can assume that the measures u^^dvolg/ fy ul^dvolg 
converges weakly to a measure ^pg(^ a{P)5p, Ci subset of Cminmin, o(P) > for all 



a{P) = 1. hm / ul^dvolg/ / u^^dvolg = a{P). 



(35) 



Hence WflCi is not empty. Say QGWflCi. Using the notation of Theorem |21 



fo = <Q) 



{2tx)- 



> 0. 



(36) 



sup Me„ 

hm inf '-^^^ > hm = /q > 0. (37) 

n->oo SUPyMe„ n->oosupyUe„ 



Lemmg^ends the proof. 



2.5 The set of limit measures 

In the previous section, we have proved that the limit measures are concentrated on the 
set Cminmin- Howcver we did not determine the concentration coefficients completely. 
Indeed let { u^Jn G N} a converging sequence with > 0, If we set = then for 

each point P G Cminmin, linin->oo t^e„(-P) = tt(-P) £ [0,1]- The concentration coefficients 
are given by 
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At least one of the a(-P) is equal to one. At present time, we do not know if the a(-P) are 
independent of the sequence { u^J\n G N}. The modulating factors a(-P), defined on the 
set Cminmin, depend only on the ratio of the relative maximum to the global maximum of 
the eigenf unction. Given a function (3 : Cminmin — >[0, 1], is it possible to find a converging 
sequence { u^^\n G N} such that Ci{P) = l3{P), for P G Cminmin? Or is the function a 
unique? So far to evaluate the concentration coefficients only the potential c was needed. 
We shall see in the rest of this section that actually only a subset of Cminmin can be 
charged. The rest of this section is devoted to study the influence of the Riemannian 
structure on this selection process. 

We will use the following deflnitions: 

Definition 4 A point P G Cminmin is called a maximally charged point if 

hmsup = 1 

(e,Q)^{o,P)maxv/„ Ue 

Definition 5 A point P G Cminmin is called a charged point if 

linisup ""'^^^ > (38) 
(e,Q)^{o,P)maxy^ Ue 

Remark. 

At this stage we do not know if indeed the analysis can be pushed further to prove that 
whether the limit measure is unique or not. For example can one flnd a double well 
asymmetric potential such that sequence of the flrst eigenfunctions will concentrate at 
only one of the two minimum points and another sequence only at the other? 



2.6 Expansion of the eigenfunction 

Recall the blown-up function: 

U!p,e{x) = = 

Ue 

By Theorem 121 (i), wp^^ converges to wp in L^(M™) when e goes to 0. 
The divided difference 

Wpi; — Wp 
U)l,e = 



as e goes to 0, converges to wi satisfying the equation: 



Lpwi + ^ CijkViUjykWp = (39) 

l<i<j<fc<r?i 
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where 



and: 



Lp = Ae 



fc=i 



Id 



c — mm 



l<i<j<k<m l<i<3<k<l<m 



k=l 



where the Cijk are constants and the Cijku C°° functions on Up. The proof of the con- 
vergence of the sequence wi^s results from various estimates that can be obtained on the 
sequence and its derivatives, which can be found in |12^ . But a rough idea of the conver- 
gence of the sequence results from the WKB expansion. The estimates are based on the 
exponential decay at infinity. Those estimates are not used in the remaing parts of the 
article. 



2.6.1 Hermite functions 



We need Hermite functions. For nG define 



where for nG 

a£ _ 2 

Kix) = e 2 — 

Note that 

Lpi7„ = 2<n, v^ACP) > 



where A(P) = (^/Ai(P), ••, VKi{P))- It is well known and easy to check that 

/ HpH,_dy = 

if p,q E 11^ and p ^ q. 

Hldy 



^Ai(P)..A„(P) 

m 

where \p\ = Y^n^^Pn, p[ = YIpJ 

n=l 

In particular:Hjjfc = -8^/X~X~>^yiyjyk, for 1 < i < j < k < m, Raj = -^Xi^/Y^yfyj + 

3 

A^XiXjyj for 1 < i < j < m , Yim = 12Ajt/j — SXfyf. To summarize, 
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Lemma 4 Convergence of the renomalization sequence. 

(i)As e goes to 0, wi^e tends to wi, solution of the equation: 



k=l 



(a) Up to a multiple ofwp 



u^iiy) + c^jkiP)y^yjykWp{y) = 

l<i<j<k<m 



H. 



(y) 



E 

l<i<j<m 



i<i<j<k<m 16^A,(P)A,(P)A,(P)< e. + ej + e,, ./X(P) > 



8AyA?(P)A|(P) 



--HeAy) + 



8^yAf(P)A,2(P) 



--HeXy) 



2.6.2 Expansion of the eigenvalue 



It is well known that the eigenvalue A^ has an asymptotic expansion of the type A^/6 + 
9e + ... (see ^301 EI])- We have already seen that: 

VUP)\P e Critic 

n=l 

In this section we want to compute 6 using minimax procedures. Recall that: 



Let us set: 



Ae = inf 

u&m{v)-{o} 



Iv 



I I T— 7 I I 2 9 

e| I Vm| + cu 



dvoln 



fy u^dvolg 



Ae — mine — JeK 
V ^ 



Then: 



e, = inf 

u€H^V)-{0} 



Iv 



dvolr. 



Jy u^dvolg 



(40) 



We shall prove that as e goes to 0, 6^ converges and compute its limit 0. It will 
follow that the limit measures are supported in a subset of Cminmin- We shall estimate 
the variational quotient 



Iv 



|2 _|_ c-mmc-y^A 2 



dvoL 



J uldvol, 
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using the previous blow up analysis. Consider a converging subsequence {u^^jn G N}, u^^ 
normalized {fy u^^dvolg = 1), such that 

uljvolg ^ 7p(5p, 

in 

where = 1 for at least one Pi. To simplifie the notations we shall drop the index n 
for En from now on. 



m 



Bp{&) 



,„ ,,9 c — mine — \feK ^ 



dvoL 



+ 



V-VJp^c^.^Bp(5) 



|V«e||^ + 



c — mine — v^A 



dvoln 



Let be a C°°function:V — s>[0, 1], 1 on V-UpgCminBp((5) such that supportcj) fl Cmin=0- 
For £ sufficiently small 



I 



,,9 c — minc — \/eK , 



V-Up6c^.^i?p(5) 

Using the relational 



dvolg < 



,,9 c — minc — \/eK 9 



(pdvoL 



,„ ,,9 c — mine — \/eA 9 



(pdvolg 



A — mmc — v^A 



A0 



u^dvoL 



Appendix 2 shows that with an appropriate constant C 

X — mine — i/eA 



V 



-A4> 



uidvolg = o(l)e™/"u 



m/4 - 2 



Hence: 



,,9 e — mine — \/eA 9 



dvolg + o(l)e'"/^M 



m/4 - 2 
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Bp{5/re) 



|V«;p,e(2/)||' + 



(c(yv^) - mmc - y/el\)wp^{y)' 



Vdetgjyjdy 



The normalization condition implies that 



uAvol 



u^dvol 



V-Upec^.^Bp(5) 
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utdvolg = o{l)e"'^W 



uldvolg = e^^^W 



Bp{5) 



Wp,e (y) y^det^(?/v^) dy 



(41) 
(42) 



Now sup \det g{yy^)\ < +00 and on each compact subset of M™, det g{yy^) con 

BpiVP/^) 

verges uniformly to 1 when e goes toO. Hence 
/igm wpiyYdy as e goes to 0. By Theorem |21 



Bp{^p/re) '^pAy)V^^^9{y\^)dy tends to 



» m 

/ wp{yrdy={2n)ff,llX~'^\P). 

^R- fc=l 

Then relations I^TE^ imply: 



in 



(43) 



fc=i 



Now we shall compute the numerator of Q{u^). 



lim QiUf) = K y lim 

e — ^ ^ e — >0 



n 



Bp (5/^) 



1 / Mi2 , {c{yy^) -mine- ^/eA)wj, {y) 



X 



detgiy^)dy.iAA) 



To simplify the notations below set 

m 

q{y) = Y.^k{P)yl-I^ 



k=l 



L 



Bp{&/m 



liv^p,e(y)||; 



{c{y^) - mine - ^h) wl,Xy) 



det g{y<fe)dy = I + II + III + IV 
+V + VI + VII 



where, denoting by V eI the euclidean gradient (^, 



1=^ / [||Vi,t/7p(y)||^„ + gw7p(y)2] -y/d^t^(^^d|/, 

VejBp(5/^) 



II 



[< Vijwp, VijWi,£ > +gwpwi_e] J det g{y^)dy 
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Ill =^ 



Bp(S/m l<^<j<k< 



^ CijkyiyjykWp\J det g{yy^)dy 



IV 



Bp (5/^) 



\\V EWi^e{y)\\Rr^ + qwi,e{yf + 2 ^ CijkyiyjykWpw 

^<:i<:j<:k<m 



he 



det g{y^)dy, 



^ =^ I XI ky iy jy kwl^ y det g{y^)dy, 



VI =/ ^ X Cijki{y^)yiyjykyiwp^e\l ^^^9{y\/^)dy, 

where ^^^(x) = 5^^' + Y:Zi=i9'''''('')''kXi.ln particular g'^'^^O) = \{Rikij{P) + RjkuiP)) ■ 
Note that 



^detg,{y) = l-^R,,{y<re) 
o 



yiyj 



where the Rij are functions on f/p and Rij{0) = RicCij{P). It is easy to see that 
limVII=/ V ^(i?ifc^,(P) + i?,fc;^(P))l/.y^^^c^|/, 



lim VI 

e — >0 



X CijkiiO)yiyjykyiwldy, 

l<i<j<k<l<m 



lim V = 

e — >0 



lim IV 

e — ►O 



|V£;U7i(y)||Km + + 2 ^ CijkyiyjykWpWi 

l<:i<:j<:k<m 

Using equation inni satisfied by wi, 



dy 



lim IV = / CijkyiyjykWpWidy 



l<i<j<fc<m 



III =^ 



Bp(V^)l<,<j<fe< 



CijkyiyjykWpdy — 



l<Z,n<m 



CijkRiniy\^)yiyjykyiynwldy. 



For reason of symmetry the first integral is and the second tends to as £ — > 0. 



lim III = 

e — >0 



43 



1=^ 



ij=l 



but 



^ / [I I VEi(;p(y) I llm + gwp(y)^] dy = [ wp < VeWp, 'n >dA 

converges to zero, due the exponential decay. Here Sp{5/^) is the sphere boundary of 
Bp{5 / \/e), 'n is the exterior euchdean unit normal to Sp{S/^) and dA the area measure. 
It is then clear that 



lim I = — 

£ >0 



^ m 

^EWp{y)\\lm + qwp{yf] - RiCij{P)yiyjdy, 



ij=l 



II = III + 112, 



III = ^ / [< VfiWp, VeWi^s > +qwpwi,e] dy, 

ye jbp{8/^) 



III = ^ / wi,e < ^EWp, n > dA, 

ye Jsp{5/^) 

lim III = 0, 

£ >0 

where an intergration by part leads to a boundary term, which converges to zero. 



< 



VEWp,VEWi,e > +qwpwi^e]YRij{y\fe)yiyjdy 



ij=l 



lim II2 = 0. 

£^0 ^ 



Finally lim Jj 



j=J\^wpAy)\\ 



2 , (c(j/ v^)-mmc— /eA)tUp_^(j/) 



^/det g{y^)dy is equal 



to lim 1+ lim IV+ lim V+ lim VI. Lengthy but straightforward computations show 

£ >0 £ »-0 £ >0 £ >0 



that with ^A(P) = ^/\JP): 



n=l 



lim I+limV = - f^- 



1 
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The computation of lim IV is more involved, where 

£ >0 



Lp = Ae + ^x^A„(P)-A(P). 



n=l 



Then 



(</)+ 



l<i<j<m 

After lengthy computations we get 



HeM + 



_4^A,^(P)A,(P) A^X,{P)\]{P) 



--HeAv) 



m 2 

J2 CijkViyjykWpw.dy - (^(^) + B{P) + C{P)) . 



l<i<j<k<m 



B{P)^ 
C{P)^ 



l<i<j<k<m 



16^\,(P)Xj(P)Xt{P) <e, + Cj + et, y^MP) > 



0<i<j<m 



E 



0<i<j<m 



m\,(P)\j{P) ^ sAKi') 



E 



Theorem 4 T/ie expansion of the eigenvalue in power of -Je is up to order three 

Ag = min c + A^/e + 6e + o{e) 



V 



where 



C'min = {P such that c{P) — mine}, 
A = inf 



J2 V^n{li)\R e C^i 

_n=l 

m 

C„,inmin = {P SUch that A(P) = yj \n{P)\P & C^in} 



n=l 



and 



— mm — , 

PeC min min ^/Xfp) 



R{P) 1 



' A.(P) 
A,(P) 



4 ^ VA.(i')A,(P) 2^ A.(P) ^ ^ ^ ^ ^ ^ 
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Remcirks. 



1. The sequence is in if" but not in if'', for s > m = dimV. This is exactly the 
critical case in the Sobolev embedding theorem. That why the limits of as e goes 
to zero are measures but not functions. 

2. The coefficients of the limit measure depends on two factors: one is the second 
derivative of the potential at some points of Cminmin and second on the coefficients 
a, that measure the ratio of the local maximum versus the global maximum. 

3. Actually, we conjecture that if two points Pi and P2 are always charged and under 
some conditions of analycity on the metric, any small neighborhood Vi of Pi is 
conjugated to a neighborhood V2 of P2 by an isometry. Such an isomctry should 
leave invariant the blow-up equation defined on the manifold and satisfies by w^. 
If in Vi, there exists a subsequence of local maximum, then such subsequence is 
also present in V2, which garantee the concentration in both places. In some sense, 
when the concentration occurs at those points for any subsequences, it implies the 
existence of hidden symmetries. 

4. Suppose that the potential has only two wells centered at Pi and P2 that are charged. 
Obviously when the coefficients of concentration are not unique, the set of possible 
limit coefficients ai,a2, where /i — aiSp^^ + ck25p^, are restricted to the intervals: 



Indeed one of the coefficient, say ai characterize the convergence of the absolute 
maximum. 

2.7 Removing the degeneracy 

In the previous analysis we have defined two local characteristics (quantities that are 
functions of g and c), Yln=i ^^niP) and —V{P) + C{P) which have properties that the 
minimum of the first one on Cmin gives the first coefficient of the developpment of in 
the power of ^/e, 



Cminmin IS the subset of Cmin where the characteristic achieves its minimum. The second 
coefficient of the developpment of A^ is equal to the minimum of the second characteristic 
on the set Cminmin- We could continue this process and define a sequence of characteristics 
X2, X3, X4, ...(Xo = miny^ c, Xi = A and X2 = 0) and subsets of Cmin: C's, ^4, C5... having 
the following property, for any integer n: 

• Xn is a function defined on C„ and depends only on the covariant derivatives of c 
and g 

• the n*'* coefficient is equal minc„ Xn- 



ai e [ci/ci + c2, 1] and a2 E [0, C2/(cl + c2)]. 




+ ... 



46 



Cn+1 = {P e Cn\Xn{P) = 



Since the set Cmin is finite, the process becomes stationnary after a finite number of 
steps i.e. C„ = C„+i and Xn is constant on C„, for n larger than n, some integer. Now 
there are two possibihties. Either Cn is reduced to a single point, in which case, this is 
the only point that is charged or the so-called degenerate case, see page 304) where 
Cn has more than one point. 

The knowledge of the expansion of Ae is not enough to compute the value of the 
concentration coefficient. It is an open problem to compute precisely this value in terms 
of the geometry and the extracted subsequence. 

Finally we conjecture that for any points Pi, P2 in Cn, there exists a isometry : fii — >■ 
f^2, where fli is an open neighborhood of Pi such that denoting by Ci, C2 the restrictions 
of c to Qi, Q2 respectively, C2 o = ci and ^^1^2 = ""^ebi ° 0- 

3 Blow up analysis with a gradient vector field 

In this section, we study the behavior of the sequence of first eigenf unctions, when the 
vector field b is the gradient of a Morse function, b = V0, with respect to the metric g. It 
is possible to describe explicitly the set of limit measures and we will see how the supports 
of these limit measures coincide with those suggested by the analysis of the Topological 
Pressure (^j). More precisely, the potential c and the field b = V0 interact to restrict 
the limits of the eigenfunctions, to certain subsets of the critical set of b. The limits 
depend on the couple of functions (c, (p). 

Somewhat surprisingly not all the attractors points are charged. It depends on the 
topological pressure (see |22]) and can be understood as followsx acts as a killing term 
(see 123113) and can destroy all the particles near an attractor; 

We begin analyzing the problem. The eigenfunction satisfies the partial differential 
equation and normalisation condition 



We denote by Sing(b) the set of singular points of the field 6. When b is the gradient of a 
function 0, these points are just the critical points of (p. 

3.1 What do we learn from the Topological Pressure 

The topological pressure can be defined as follows: 



(46) 




(47) 



m 



Pr 



inf 

P&Sing{b) 




(48) 



i=l 
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where Aj(P), 1 < i < m, are the eigenvalues of the hnear part of the field gradcf) at the 
point P. The definition for Morse-Smale fields can be found in (j2D]); for more general 
vector fields in (|22|). Because the infimum in ()48|) is taken over a finite number of 
points, it is attained somewhere but not necessarily at an attractor (that is clear from the 
formula) . 

It will be useful to perform a gauge transformation: 

b = V(j). 

By this gauge transformation, equation ()46|) is transformed into 

2 4e 

with the new condition of normalization Jy v^dvolg = 1. We will analyze the set of limits 
of the measures v^dvolg as e tends to 0. 



3.2 blow-up Analysis 

It has been proved in the second section that the sequence sup^ converges to infinity as 
e goes to zero and the sequence converges uniformly to zero on every compact set that 
does not intersect the singular set of the field. As in the field- free case, if Pg Sing(6) we 
choose a normal coordinate system (xi, ...x^):!! — >R, centered at P, defined on a domain 
U such that: 

1) xix. ..xxm{U) contains the closed ball Bp{6) centered at P and having radius 6 > 0. 

2) for all i, j, 1< i, j< m, agr(P) = A,(P)5,,. 

3) UnSing(b)={P} 

On the magnified set -^xiX...xx.m{U) , we can define the function We:We{y) = ""^^^^^ 
where W£=max Vs- 



The main theorem of this section is the following: 

Theorem 5 • (i) Suppose that the vector field b is gradient-like, b = V0, where has 
the expansion near each critical point: 0(P) = 0(-P) + Yl^i Aj(-P)(a;')^ + 
The set Ai of all possible limit measures of the sequence v'^dvolg, is described by: 

V = {fie M{V)\i2 = cp6p, cp > 0, ^ cp = 1} 

Pes Pes 

S is set the critical points of b, where the Topological Pressure is achieved. 

• (a) When all the points of S are maximally charged, the coefficients cp can be 
computed explicitly, cp = ^ ^^rl^' i wpm-i/2 
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When a partial set S of S is maximally charged, the formula changes to: 



nr 




X]pe5 111 


|Ai(P)|-V2 



Cp 

Lemma 5 For each e, let us denote by A^^, the set of all maximum points of in the 
manifold. There exists a constant A depending only on c and b, such that: 

sup c?g(g, Singib)) < Ae 

It follows from this that the set of limit points of the set of maximum points of is 
contained in Singib). 

Proof. 

If we denote by a sequence of maximum points of v^. Using the Maximal Principle, 
applied to equation pHjl . it is possible to estimate the velocity at which the sequence of 
maximal points converges to a critical point of the field. Indeed, we obtain that 

(e(Q.) + ^+l™%A.<C 

for some constant C. This implies that | |V0((5e)| p < 4Ce. Now it is easy to see that there 
exists a constant F depending only on c such that for P G V, dg(P, Sing(b))< r(c(P)-min 

c). Take A=4rA. Then dg{Qe, Sing{b)) < Ae. The lemma shows that Qe tends to Sing{b) 
at least as fast as ^/e. In fact, theorem below shows that it tends faster. 

Theorem 6 Concentration with a gradient vector field. 

(i)For any Pe Cmin,cmy sequence of Vs's, with e tending to 0, contains a subsequence 
{vs„}, such that the sequence of blown-up functions _ at P converges to a function 
w;IR™'->]R_i_, both in the L'^norm and the CP° topology. 

(a) w satisfies the equation: 



w = Xw (49) 



< w < maxw = 1. 



(iii)When the eigenvalue A is equal to the Topological Pressure at P, the solution w is 
nonzero and given explicitly by: 



n 

in particular, 



w[x] = 

i=l 



m 

|A,(P)|x2/2 



w^(x) 
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Moreover, 

lim e""/^ sup vl = K{c, 0) 



where K{c, (f)) is a positive constant depending only on the functions c, 0. 

(iv)If all the singular points at which the topological pressure is attained (set denoted 
by "top") are maximally charged, 



E{Ki|A.(P)|-^/^|Petop} 
Moreover 

€^0 e^i^ 



Remark. 

At this stage, let us point out the difference with the pure potential case, studied in the 

last section. There the scaling factor was and the concentration took place on certain 
minimum points of the potential only. In the present situation, the scaling factor is \/e 
and the concentration is determined by the couple field-potential ( V0, c). 

Before proving the theorem, two lemmas are needed. The first one provides estimates 
of the first eigenvalue Ag and the second, estimates of the decay of renormalized sequence 
w^{x). 

Lemma 6 The first eigenvalue satisfies the inequality 

Q<K<min{c{P)- ^ \i{P)\P e^\ng{h)} 

{i,Ai(P)<0} 

where Sing(h) is the set of critical points of the field b — V0 and Xi{P), I < i < m, are 
the eigenvalues of the field at the critical point P. 



Proof. 

This is a consequence of the variational approach. Similarly to the case without field, 
define the variational quotient 



JyV^dvolg 



where c, = c + + Then 



A, = inf QJv) > 

veH^{V)-{Q} 
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In the neighborhood of a singular point P, we have an expansion, 

c,(x) = c(P) + ^E^^ + J2 a,{P)x' + + c,jk{x)xiXjXk 

i=l i=l <ii<j<k<m 

where the Cijk are C°° functions on Up. In order to get an upper bound for the eigenvalue 
Xe, we estimate the quotient Qe, for the following test function with compact support, 

= 0, otherwise 

where /ij = 1^^^^ and A/p(p) is the connected component of the set {x|^"^]^ < p} 

containing and p is taken so small that Mp{p)n {x|^™^a;f = 6^} is empty. So A/p(p) 
is contained in Bp (5). After computations similar to the ones in Lemma |21 of the section 
121 we obtain 

Q,{^,) = c{p) + + 9 + -^h^ + ^^^) 



i=l 



i=l 

Using the fact that '^^^^^^'^ = — Y^^=i ^i{P)i we obtain that 

Q.(^.)<c(P)+ Y ^«(^)- 

X^{P)<0 

Since a similar test function can be built in the neighborhood of every critical points, we 
obtain the estimate of the lemma, 

A, < min ^^Jc(P) - Z{^^(^)|1 <^<m, A.(P) < 0}]. 

{P&Sing(P)} ^ — ^ 



Lemma 7 The blown-up sequence w^, defined by 

wAx) = 

satisfies the following properties: 
For all Eo g]0, 1[: 

1. sup Jbp{s/ i^)'^£iy)'^^y < +00. More generally, for any continuous function f:[0, 

]0,eo] 

IJxMJ^ — >M.,{e,y) — >/ (e, y), having at most polynomial growth at infinity in y 
uniformly m e, sup jBp{s/^) fi.^^y)'^e{yfdy < +oo. 

]0,eo] 
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2. the set of restrictions w^\Bp{S/ y/e) , e g]0,1], of the to the balls Bp{5/^) sat- 
isfies the following condition: for any t] > 0, there exists a compact Kc and a 
e{ri) > such that 



Bp{5/y^)-K 



for all e g]0, e{7])\. 



Proof of Lemme(71 

The proof is similar as the one given to prove the decay estimate in the case of the pure 
potential case and is based on the Feynman-Kac integral representation of the solution. 
The Feynman-Kac formula gives that 

where we recall that q = c + + ^^^^ and is defined in relation (j^ .Tf has the same 
meaning as in Section 12)). Since the function c+ is bounded, it is enough to estimate 

Ex{e~ -^0 Using the scale change x — >■ e^/^x, this term can be estimated as in 

Proposition 1211 and the results are similar. Indeed, because ^^^^ = J2^=i ^niPYUn-, 

the estimate is valid in the neighborhood of any critical point. | 

Proof of the theorem. Consider the renormalized sequence w^^ 

{y) = satisfies 

the partial differential equation: 

Ag^w,{yVe) + \c{yVe) + + ^ 1 w,{yVe) = Kw, (yVe) 

-^Xi X ... X Xm{Up) (50) 
< We < 1 

To study the equation EOl note that the coefficients of the partial differential equation EOl 
converges uniformly on every compact set as e goes to zero: ge{x) = g{yy^) converges 
to the Euclidean metric (see the previous section for more details about the type of 
convergence, usually in C^'" of any compact). 

Since the coefficients of the partial differential equation EOl remain bounded as e goes 
to zero. Indeed, due the velocity of convergence of the sequence QJHl 

n=l 



m 
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Classical elliptic estimates show that converges uniformly on every compact set of M"* 
to a solution w of 

AEw{y) + (^c{P) + ^(P) + f2\l{P)y'}j w{y) = \w{y), on (51) 
< w < maxw < 1. (52) 

Let Q* be a limit point of the sequence since converges uniformly on every compact 
to w. We(^) converges to w{Q*). Hence w{Q*) = 1. Because w belongs to L2(M") by 
Lemma[7|and is not zero, w an eigenfunction of Lq in L2(M"') associated to the eigenvalue 
A. 

Moreover, using the results of lemma 

A-(c(P) + ^(P))<f;|A„(P)| (53) 

^ n=l 

We will now prove, in order for the function w to be nonzero, necessarily A is equal to the 
Topological Pressure. 

In order to prove that /i = A — (c(P) + ^(P)) = Yl'^=i l'^"(-P)|, one must study 
the spectrum of the self-adjoint operator defined by Lp = A^; + J2n=i ^"^^^^" in L2(M"'). 
This operator is compact and has a discrete spectrum. For all these, see ^7\. Using the 
arguments in the proofs of the main theorem in Section |231 we know that the eigenvalues 
and eigenfunctions of the operator defined by Lp in L2(M"') are respectively 

X:(2a. + l)^4^,a6Z'^ 



and 



Ha{x) = llh^Me~^^''^' where h,ix) = l'^^{e-^' 1^}. 



i 

The lowest eigenvalue of Lp is 

n 

1 

Inequality ESI shows that the eigenvalue A — (c(P) + ^(P)) of Lp in L2(M") is at most 
equals to Ai. Hence it is equal to Ai and the corresponding eigenfunction it; is a positive 
multiple of Hq. Since the only maximum point of is 0, = and w = Ho{0). 
In particular, the previous analysis shows that 

limAe = (c(P) H —(P)) + Ai = Topological Pressure. 

£—►0 2 

This implies that a point P is not charged if the Topological Pressure is not attained at 
P. 
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Using the exponential decay of the last lemma, at the point P, where the Topological 
Pressure is attained, we have 

vidvolg = e'^/'^vl / wldvolg^. 

Since the sequence converges in L2, this implies that the sequence e^^^vl converges 
and we get, 

lime'^/V = C(c, 0). 
The concentration coefficient can be computed, if we denote 

cp = lim / v'^dvolg = K{c, (f)) / Wq. 

Using the L2 normalization conditions, the Blow up analysis and the fact that the sequence 
converges to zero on any compact sets that does not intersect the critical points where 
the topological pressure is attained, we get. 



l = ir(c,0)5^ / I 



wl + o{l) 



which proves the last part of the theorem. We shall omit the proof of statement (iv) 
which is similar to the one given at the end of sectior|21 

Remarks: 

1. Can one find explicitly the asymptotic expansion of A^, 



A. = J]c,e^/'. (54) 



fc=0 

2. We proved earlier that cq is equal to the Topological Pressure and captures enough 
information to locate the supports of the possible limit measures. This is surprising 
because it was not true in the pure potential case. See the section (j2.7p 

4 The general case for the first eigenfunction prob- 
lem 

In this section, we shall consider a special class of vector fields h such that the limits of 
the first eigenfunctions concentrate on the limit sets of h. In particular if the limit sets 
of the vector field contains limit cycles, we prove that the weak limits of as e tends to 
zero, concentrate on these limit cycles, when a 

Here we will consider the following class of vector fields : h = gradgC + Q where is a 
Morse-Smale vector field and £ is a special type of Lyapunov function of Q. We will study 
the behavior of the first eigenfunction for the operator the drift of which is belongs to 
this class. We will need the following lemma proved in the appendix. 
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Lemma 8 Given a Morse Smale field Vt there exists C°° Lyapunov functions for Vt, taking 



the value zero on all the repulsors ofVt and such that the function = idl^'^llo + i < 



grade, VL >g) is positive except on the recurrent set of Vt where it is zero. 

Definition 6 A Lyapunov function with the property stated in lemma\^will be called a 
special Lyapunov function. 

The main theorem can now be stated: 

Theorem 7 On a compact orientable Riemannian manifold Vm, consider a Morse-Smale 
vector field Vt whose recurrent set consists of the stationary points Pi,...,Pm and of the 
periodic orbits Ti, .., Vj^. C denotes a special Lyapunov function. 

For e > let X^, denotes respectively the first eigenvalue and the associated eigen- 
function of the operator 

= eAg + W + c , onVm (55) 
Then the weak limits of the family of probability measures 

e-'^/'uldVg 



are of the form 

N M 

j=i i=i 

where the measures fij are supported by the limit cycle Tj and the Pi's 1 < i < N are the 
critical points ofb. the 7p. are scalars. 

N „ M 



7 = 1 -^^j i=l 



Remark. 

In a future paper fT?|, we shall prove that for an appropriate choice of Lyapunov 
function, no concentration can occur on a set of isolated points on the cycle, 
point where the topological pressure is found 

4.0.1 Proof of theorem 2 

The proof of theorem 2 will occupy the next two sections. Consider the following trans- 
formation u — > V. : 

V = e~^u (56) 
u = e^f (57) 
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and the second order differential operator L'^: 



It is easy to see that: 



AgU = e2<= 



(A,„-1<V^.V.„>,+JM-^I» 



and 



Finally, 



< Vm, b >g= 



< \/v,b >g +v 



< V(f),b>, 
Ye 



L'M = e* 



eAgV+ <b - V0, Vv >g +v \ c + 



<V<f>,b>g , Ag<f) IIV0II 



2e 



4e 



Setting b^ = b — Vcp, we obtain : 



(58) 



eLv = e^AaV + e<bs, Vv >„ +v sic + ^ + ^ ^ 



IIV0II: 



If we set Ve = e~2iu£, the equation 

eAgUs+ < b, VUe > +CUe = XeUe 

becomes : 

e^AgVe + e < 6<^, Vf^ >g +CeVe = eXeVe, on V (59) 

where c, = e(c+ f ) + + 

We set = fi. c — ^ — ^^^^ > 0. Taking = a special Lyapunov function for Q, 
in the previous paragraph, we will prove that eA^ tends to the minimum of the function 
"^c = ^^f^ + when e goes to zero. 

Lemma 9 Under the assumptions of Theorem^ on the vector field 

lim eA, = = min\E'/: (60) 



Proof: Multiplying equation ()59|) by and integrating on V, we obtain: 

edw{b^) 



[e'\\^v,\\t + 



+ c^)vf\dvolg = eAe / vldvolg 
Jv 



(61) 
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Normalize so that Jyv'^dvolg=l. From equation (jHH), we obtain that hme^^oe-^e > 
min^^. An upper bound of the function eA^, in the case when = £, will be computed 

using the following two results: 

1) Let 5 be a regular sub-domain of V and X[{S) be the first eigenvalue of L'^ on 5* for 
the homogeneous Dirichlet problem. If there exists a positive number A and a positive 
function ip on S, zero on the boundary of S, such that L[ip < eAip, then X'^{S) < eA (see 

my 

2) If Ae is the first eigenvalue of Le on V, then eX^ < X'^{S).l^ote that eX^ is the first 
eigenvalue of on V ([3^). 

We proceed here with the proof in the case of a limit cycle and leave the case of a 
singular point, the easier one, to the reader. We construct a test function as follow. Let 
be a stable limit cycle of b. Let be the connected component containing uj, of the 
set {P G V I C{P) < C-iuj) + 5}. For 5 sufficiently small, is a regular sub-domain of V 
and a neighborhood of uj not intersecting any other recurrent set of he- For simplicity we 
set he = B. 

Define the function : F^— > M, as follow : 

Denoting by V the gradient with respect to g, we have: 



-C/{2e) 



eA,C 



eA,C 



2 4 



+ 



iv/:i 



Let i? be a positive number. 



L' jp — eRip = e 



-5/(2.) 



e c 



R){e 



S-C 
2e 



|V£| 
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Below we shall determine R and 6 so that L'^ip — eRip < 0. Assume this done, using results 
1 and 2 above, with S' = F^ we get that eA^ < A'^(S') < eR. Hence lim^^o ^X^ = 0. We use 
Lemma ITT] to determine R and 6. We get that up to terms of order four at least, using 
appendix I, 



40 - ei?0 = e-*/('^)(e(c - R){e^ - 1) - etVgA - (1 - /i/2) \\Ax\\') 

We set 6 = 2ae where a will be determined later. There are two cases. If £ < ae, then 
— 1 > — 1 and we can choose a so that (e? — 1) is arbitrarily large. < Ce 

where C > depends only on the matrix A and a. When e converges to zero, for R big 
enough and C < ae: 



L'^(f) - eR(j) <e[{c- i?)(et - 1) + \trgA\] < 
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When £ > ae, then 

(c — R){e^^ — 1) < for R large enough 

and: 

L'^(j) - eRcj) < e{trgA - (1 - ^)aia). 

Because C{x) > ae, it imphes that ||Aa;||^ > aiae, where ai depends only on the matrix 
A. Finally L'^cj) - eRcp < if < a. Hence for R large enough, L'^cj) — eRcp < in 

r^. Using the results of (0|), we obtain that eX^ < XeiXe) < ^R- This ends the proof of 
Lemma El 

Now we will prove that any weak limit of v'^dvolg/ jy v^dvolg as e goes to zero, con- 
centrates on the limit sets of h. 

Lemma 10 All weak limits of measures v^dvolg/ jy vidvolg as e goes to zero, are con- 
centrated on the minimum set of the function \E'£. 



Remark. 

By construction ^ c. attains its minimum on the limit sets of VL only and this minimum is 
zero (note that on a repulsive set, L is maximal, and ^ c is also zero). 

Proof of the Lemma IIOL 

Proceeding as in the self-adjoint case, we multiply equation ()59|) by f ^0, where is an 
arbitrary test function. We obtain after integration by parts: 



V 



'v 



Ce converges to and eA^ to = niin\E'£ as e goes to zero. 
Let /i be a weak limit of the measures v'^dvoL then 



j (f){imn^c - ^c)d^i = 



(62) 



Since equation (jU^ is true for any test function cj), dfi = on the open set where 
Hiin\E'£. This shows that the support of /i is contained in the limit set of Q. Let 5*1, ...Sp be 

the stationary points of Q and let Fi, ...r„ be the limit cycles. We have the decomposition: 



k=l 



where /i^ is a measure supported by F^,. 
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4.0.2 Remarks on the limit measures 



We end this section with two remarks. An interesting question is whether the repulsive 
sets can be charged by limit measures. Were it to hapen this is similar to the fact, well- 
known in the large deviations theory, that the trajectories of a stochastic system can 
exit with a positive probablity the basin of attraction of an attractor of its drift. The 
boundary of the basin plays the role of the repulsive sets here. In other words a particle 
moving according to such a system can escape from the basin of attraction. 

The second remark deals with the relations between the limit measures whose existence 
was shown above, and the measures, called equilibrium measures, which maximize the 
topological pressure P, defined in the introduction. More precisely P is given by the 
formula: 



sup 



V 



(63) 



where M is the set of probablity measures invariant by the flow of the drift. A measure /z G 
M is an equilibrium measure if it maximizes hy{F^) + Jy{c+(f)'^)di>. For more explanations 
see(|22])- 



4.1 Final Remark 

The characterization of the set of possible limit measure in the nonvariational case is 
studied in jT3] , where we prove that the concentration can occur along some submanifolds 
of the recurrent sets of a hyperbolic field b. In some cases, we are able to show that 
the measure is absolutely continuous with respect to the Hausdorff measure induced on 
the set. The order at which the Lyapunov function vanishes, in the neighborhood of the 
recurrents, play the role of a filter that allows the eigensequence to concentrate only along 
a subset where the topological pressure is achieved. 



5 Appendix 1 

In this appendix, we give the construction of Lyapunov functions for Morse-Smale vector 
field satisfying the condition i( ||Vi2||^ + 2 < VC,b >g) > 0. We start by a local 
construction near the recurrent sets and then give a global construction on a compact 
Riemannian manifold. 



5.1 Local construction of Lyapunov functions 

Lemma 11 Given a Riemannian manifold (V, g), of dimension n and a vector field b on 
V, for any hyperbolic stationary point of b or hyperbolic periodic orbit there exists a local 
Lyapunov function C at that point or periodic orbit such that on the domain of definition 
of C outside the set of hyperbolic points and periodic orbits, "^{L) = i( ||V£||^ + 2 < 
VC,b>g) > 0. 
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Since is a Morse-Smale vector field, the limit sets are contained in the union of critical 
points and limit cycles. We refer to Kamin (see ^3 for the construction of L in the 
case of a critical attractive point of the field. 

First we construct L near an attractive orientable periodic orbit F. Then the normal 
bundle A/" to F is trivial. Denote by T the minimal period of F. Let p:]R ^ be a 
periodic trajectory of period T, the image of which is F. We consider a Fermi coordinate 
system (6', Xi, .., in a neighborhood U. of the cycle, d is the cyclic coordinate and 

(ei(^), .., e,i_i(^)) is the corresponding orthonormal frame field trivializing the normal 
bundle of F i.e.dx*(ej) = the Kronecker symbol, 1< i,j < n — 1. If ^ is a point in the 
domain U of this coordinate system, 

n~l 

^ = exPp(e(0)Ea;i(Oei(^(0)] 

i=l 

In these coordinates: 

n-l 

g = gijdx^dx-^ , 

i,j=Q 

where we set dxo = d9. Along the cycle, gij{9,0) = 6ij, I < i,j < n — l,goj(6', 0) = 
gjo{0,0) = 0,1 < j < n — 1. The Christoffel symbols associated to the coordinates 
(xi, .., a;„_i)are zero along F : T^j{9,0) = 0, for i,j,k G {l..n — 1}. In a neighborhood of 
F, gij{6, x) = 6ij + 0{d{x)). In this coordinates the equation ^ = fi(^) can be written as 
follows: 

^ = 1 + 0{d{x)) (64) 
X = B{e)x + 0{d\x)) (65) 

The solutions are given by 

e{t) = ei^t)), x,{t) = xMt)), i<i<n-i, 

where is trajectory of b and d(x) is the distance from x to F, B{6) is a (n-l)x(n-l) 
matrix-valued function. 

Consider the solution X:M GL[n — 1;]R] of the matrix equation X = B{6)X such 
that X{0) = Idn~i- Then there exist a (in general complex) matrix D and a matrix 
function P:^ G M ^ P(^)eM[(n-l)x(n-l);C], such that X{e) = P{9)e'^^ for all 9 e R. 
Because the orbit F is hyperbolic and attractive, the real parts of the eigenvalues of D 
are negative: the eigenvalues of D are exactly the characteristic multipliers of the orbit. 

Denote by D* the complex conjugate transpose of D. Let /x be a strictly positive 
parameter. Consider the matrix Lyapunov equation in the unknown matrix A: 

AD + D*A = -/iA^. 

It is well known that there exists a symmetric positive definite (n-l)x(n-l) matrix A 
solution of this equation, given by the formula: 

r+oo 

Jo 
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We define a local Lyapunov function — > IR as follows: if ^ eU; 



n-l 



E^^(Oe.(^(0)], (66) 



k=l 



then: C — (AX , X)^(n-i)wh.eice X =P(^(^) ^(xi(^), ....,X(„_i)(^)) and the bar denotes the 
complex conjugate. Let tG M+ .^(t)eV, be a positive semi- trajectory of b contained 
in U and let X(t)=P(^(0)-^(xi(e(t)), X(„_i)(e(t))). 



dcm) 



dt 

dmt)) 

dt 



{AD + D*A)Xit),Xit)) +0{\\X{t)\\') 



-X{AX{t),AXit)) + 0{\\X{tm 



where for U, Ve C^" , (U, y)=J2k'=i^ U^Vk. Hence L is decreasing along the trajec- 
tories of the field b in a sufficiently small neighborhood of V. 



(ycm).mt))) = {Ax{t),Ax{t))+o{\\xm 



{vcmi'^mt))-^{mtimt))+o{\\x{tm 

and 

A£ = -2ir(A). 
Recall that *(L) = ^ + and thus 



^>{) = {l-Xl2){AX{tlAX{t)) + 0{\\X{ 
For 1 — A/2 > 0, '^{L) is strictly positive in the neighborhood of the cycle. 

Assume now that we have a periodic orbit 7 of b with minimal period T which reverses 

the orientation. Let 11 : V" — V be the covering of V associated to the cyclic subgroup of 
7ri(V) generated by 2(7] ([7] = homotopy class of 7). This is a Galois covering space with 
group Z2. Denote by i the nontrivial deck transformation of this covering, g and b have 

unique hfting to V still denoted by g and b. r=n~^(7) is a periodic orbit of the lifting 
of b (i.e) b of minimal period 2T. We can apply the preceding construction to F but with 
extra care here because of the deck involution i. g and b are invariant by i. We choose 
the neighborhood U invariant by i and the Fermi coordinate system so that: 

Tiek{e) = ek{e + T) for < k < (n - 1) 
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Then the previous construction produces a Lyapunov function L for the hfting of b 
in a neighborhood of F.This function is invariant by the Galois group of the covering and 
hence can be pushed down by H to a Lyapunov function L of b in a neighborhood of 

7:L=L o n. 

For repulsive orbits the same construction applies changing b into -b. Finally the 
general case of a general hyperbolic periodic orbit F can be easily handled by noticing 
that F has a basis of open neighborhoods U such that U is diffeomorphic to the fiber prod- 
uct over F of ZYnW*(F) and WnW"(F) and then "patching up" the Lyapunov functions 
constructed above for the restrictions of b to the stable and unstable manifolds W*(F) , 
W"(F) of F. 

More precisely, due to the transversal intersection of the stable manifold with the 
unstable manifold at the hyperbolic set, a Lyapunov function is built in each set as follows: 

(resp. iV" ) denotes the normal bundle of the stable (resp. unstable) manifold at F. 
Inside each normal fiber, a Lyapunov function is found. The coordinate system at a point 
^ eU is such that ^ = eWp(ei^)[Y.kZl •«fc(0efc(6'(0)] = ^Wp{9{ol^s + then as in the 
first part of the lemma, we can defined find two matrices Au and Ag solving the matrix 
Lyapunov equation in each subspace A^" and A^* respectively. If we denote n„ = dimN^ 
and. Us = dimN^, then Uu + Us — 1 — n and we can now define the Lyapunov function C 
by : = (AsYs, y;)R(n.) - (AyYy,, y„)]R(n„) where Y„ = P„(6l(^)"^(xi(0, ....,X(„„)(0) and 
Ys = Ps{9{$,)~^(xn^+i{^), ....,X(„_i)(^)). Here P„ and Ps are two matrices defined as in the 
first part of the Lemma for the solution of the matrix equation in each normal bundle. 
From the construction, it follows that C is positive and decays along the trajectory of the 
vector field. 

We can repeat this type of construction locally in a any neighborhood of a limit set 
of the field. Indeed since the limit set is hyperbolic, the construction of the function jC 
is possible on each attractive or repulsive fiber and due to the transversal intersection we 
can match the construction and the function £ vanishes only on the limit set in these 
neighborhoods. Outside the hmit set, there is a large choice of extensions of the function 
C. 

5.2 Global construction of Lyapunov functions 

In this section we prove the existence of a global Lyapunov function, assuming the exis- 
tence of a local one (proved in the last section) . 

M is a C°° compact Riemannian manifold and (pt, a Morse-Smale fiow. F is the vector 
field on V, generator of the fiow. Let us denote by Q the limit set of the fiow. Q is endowed 
with a partial order as follows: if tui, a;2belong to Q,u!i >- uj2 if there exists a trajectory 
the a— limit set of which is cui and the a;— limit set is 002- This partial order determines a 
filtration of Q: ft — D,o D fti D ... D fl^ as follows: Let fto^max be the set of all maximal 
elements of Q. Set = Q — ilo,maa;- Let fii^max be the set of the maximal elements of 
Qi. Set Q2 = ^^1 — ^^i,max and so on. When Qn has been defined, let ^ln,max be the set of 
all maximal elements of fin- Set fln+i = ^n,max- This filtration ends at a certain m. 
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Qjn is the set of all minimal elements of Q. 

5.2.1 Local Lyapunov functions 

A strict Lyapunov function for F on an open subset O of M is a C°° function C : O ^ 
M such that dC{F) < in O — . For uo E Vt there exists two relatively compact 
neighborhoods O^; and of u in M„ such that the closure Ot^oi in M, is contained 
in and a strict Lyapunov function C^^ : O^^ — > M such that: 

(i) n d^' = 0, for a;, a;' e Q, a; 7^ u' . 

(ii) The boundary (90^^. is the union of dsOi^^UduOt^UdiOt^ where the three components 
are compact codimension 1 submanifolds with boundaries ddgO^, dduO^^, ddiO^ such that 
ddsO^ U dd^O^ = ddiO^ and dsO^ n d^O^ = 0. 

(iii) dgOij and duO^ are transversal to the flow and diO^is foliated by arcs of trajectories 
of the flow (f)t linking dgO^ to duO^j. 

(iv) is constant on dgO^^ and d,jD^. 

(v) dsO^ n W^iuo) (resp. duOu) H W^iu) ) is a compact sub-manifold contained 
in l^*(cj) (resp. W^{uj)) which is a section of the restriction of the flow to W'^{uj) (resp. 
W'^{ijj)). If u) is maximal in VL, diOi_o = = c^sOt^. If uj is minimal in fl, diO^j = = duO^j. 

5.2.2 Construction of the global Lyapunov function 

The construction is inductive. To fix the ideas, we can always assume that, if is neither 
maximal nor minimal, jC^{dsOi_^) — 1 , jC^{duO^) — —1 jCi_j{u!) — and jC^{Oi_^) — [—1, 1]. 

In the case u is maximal , we can assume that O^j is an open ball whose boundary 
duOaj{C^{uj) = } is a sphere that Cuji^^) = —1 and Cuj{Ouj) = [—1,0]. 

In the case u is minimal , we can assume that O^^ is an open ball whose boundary 
dsO^^ is a sphere {C^{u) = }, that ^^^(a;) = 1 and >C(^(Ou;) = [0, !]• 

The construction starts as follows: Let Mq — {\jO^\u; e Jlg^"^}. Mq is a compact 
manifold with boundary{ \JduOoj\uj G 0,'^°'^}. Define jCq : Mq R, as follows Co\o^ — JC^,- 
It is clear that 

(i) Monu{O^,cjGQi} = 

(ii) Mq is a sub- manifold of codimension with boundary dM^. 

(iii) >Co ■ is a strict Lyapunov function on Mq and £0(^-^0) — ~^ 

Assume now that we have constructed a sub-manifold M„of codimension in M with 
boundary (9M„, and a strict Lyapunov function : M„ M, such that 

(i) Mn contains U^^gr^^^'^ in its interior and M - M„ D U{Ouj\(jJ G 

(ii) Ln is a constant on c?M„ and equal to c„, say. 

For each u e 1]™^^ let K = dM^ D U{0_t(a,O^)|t > 0}. Then is a compact 
sub- manifold of dMn of codimension with boundary — dMn n {(j)^t{ddsOij)\t > 0} 
in dMn- 

Also 

(i) n V;,/ = for cu, J e and u ^ J 

(ii) K) contains the manifold dM^ U in its interior. 
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There is for each lo G f^^+i a function T^^, : 14 ^]0,+oo[ such that for any 
X e K;, ^ 0^,0 < t < T^.,i^(a;) and 0(T^,^(a;), x) G a^Of^. 

Let e be a smaU positive number e < infx(zv^Ts^uj{x) and let C M x be the subset 
{{t,x)\x e K), —e<t< Ts^uj{x)}. 4> maps diffeomorphicaUy into M„. and (p{Dlj) 
are manifolds with corners. 

We can choose a C°° function : D]^ — >]0, +00 [ such that 

(i) if X e and -e < t < 0, c^{t,x) = -iCn{4>t) 

(ii) if x G K) and Ts^^{x) - e <t < Tg^u^ix), c^{t,x) = -f^C^{(t)t). C^^ is defined on a 
neighborhood of dgOa- 

Let 

T— sup / Cu,{t,x)dt. (67) 

Define now a function A^^ : K; — >-]0, +00 [ by the formula Xoj{x) = 'rlZ^CV" dt It is a 

r-Jo"'" cu,{t,x) 

C°° function. 

Let u be a C°° function R — i^jO, +oo[ such that suppw g]0, 1[ and u = 1. Let us 
define the function a^^ : +oo[) as follows a^(t,x) = Xi_^{x)u{j; — ^p^)- ^ 

function. Define c^^ : D^j -^]0, +00 [ as the sum c^, + a^, then 

(i) /o^--^")c^(t,a;)dt = r 

(ii) c^{t,x) = —-^Cn{(t)t{x)) for all x G K; and — e < t < 0. 

(iii) Coj{t,x) = -f^C^{(j)t{x)) for all a; G K, and Ts^^{x) - e <t < Tg^ui^x). 
(iii) imphes that for all x G K; all t such that Ts^uj{x) — e <t < Ts^u>{x), 

/ c^{s,x)ds^ / c^{t,x)dt -\- / c^{t,x)dt - i —C^((l)t(x)) (68) 

-'0 Jo Jts,^{x) Jts,^(x) dt 

= + C^{(f)Ts,u,{x),x) - C^{(t>t{x)-,x) = + 1 - C^{(t)t{x),x) (69) 

For each uj G Vf^^l let = (i>{Dl) U O^. Define a function : ^ R as follows: 
if y G (f){Dl),y = 0t(x),-e < i < 7;,u;(a^), J^M = c„ - /J Ca;(s, x)ds. If y G O^^, 
J^ojiu) = c„ — Fj^ — 1 + Coj{y) It is easy to see that jC^, is a strict Lyapunov function on Uui 
and A;(K,) = c„. 

We will now proceed with the construction of M. For each lo G f^^+J, let CV^ be a 
collar for dV^ in K; such that W'{V^) n aM„ C - CV^. There exists a function 
Tu,w ■ CV^ -^]0, +oo[ such that (f){t, x) ^ dJD^ if -e < i < Ty_^^ and (^{Ty_^^{x),x) G (?„0j^. 
Let 

iv„ = {dM,, u{VL,cc. G 1]™"?}) u {cv^\w e nZi) 

Nn is a submanifold of codimension with boundary in dMn. It is easy to construct a C°° 
function : A^^^ ^]0, +oo[ such that T„ = T„,^ for any u; G O^^^J. Let i:)" = {(i, x), x G 
NnO <t<Tu{x)}. (f) maps diffeomorphicaUy on a manifold with corners of dimension 
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in M. Then we can construct a C°° function c : — >]0, +oo[ such that: 



c(t, x) — 



dt 

Let 7 : iV„ — >]0, +oo[ be the function 7(0;) = /q^"^'^'* c{t,x)dt. 7 is a C°° function. We 
define c : — ^-jO, +oo[ as follows: 

e(,x)^£||(r.2, 
then c is a C°° function and /^f"^''^ c(t, = T + 2. If a; e Q^^^f and x e CV^, 

7W = / c{t,x)= I ^ £^{x) - £^{(j)T^{x)) (70) 

7(2:) = c„ - (c„ - 1 + r + C^{(l)tTM) = r + 2 (71) 

dC^{(j>t (a:)) 







Hence for a; e VL'^^l, c e CK;, -e < t < T^{x), c{t, x) = c(i, x) - . 
Let 

M„+i = M„ U 0(1^") U {U[/^|o; e n^^}. 

Define : Mn+i ^ M as follows: if |/ G M„, £n+i(|/) = jCn{y) if 1/ ^^-D", y = 4>t{x),x e 
N^, Cn+i{y) =Cn- jlc{T,x)dT and \iy eU^, u: e ^n+i-, ^n+i{y) = 

This function is well defined and C°°. If y G (t>{D'^) n^^^ornax LiU^, then y — (j)t{x),x e 
CK, for < i < 7;(x) = T^,u,{x), but then 



Cn - y c(r, a;)(iT = Cn + J ^^'^^^^^^'>^ = c„ + £^(0t(a;)) - £^(a;) = £^(|/) 

By construction, Cn+i is a strict Lyapunov function on M„+i. M„+i is a manifold 

of codimcnsion in M with boundary 9M„ = U{duO^\uj E ^^^"i} U N where N = 
{(f)(Tu{x) , x) , X G Nn}- If y G duOi_j nN then y = (f)(Tu{x) , x) , for some x G CK;.ru(a;) = 
Tu,u>{^)- We have 9M„+i = C~l^{cn+i), where c„+i = c„+i - F - 2. | 

Finally, we remark that b MS is equivalent to O, MS, when 

b = n + vc. 



5.3 Appendix 2 

In this appendix we are going to provide the estimates needed in (*). On the Riemannian 
manifold (V, g) consider an operator of the form: 



L^sAg + ee{h) + (V' + Vec^), 
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where ip is a non negative function on V and c^, e G [0, 1], is an arbitrary function on V, 
parametrized by e, smooth in all the variables . 

We assume that we have a familyjuele g]0, 1]} of smooth functions on V such that: 
Ue > for all e and 

eAgUe + e9{h)ue + (V" + \/^Ce)Ue = 
Then we can state a lemma: 

Lemma 12 There exists constants C, eQ>0, depending only on g, h, ip, and max |ce| 

[0,1] xV 

{but independent ofe\) and a universal mapping 'y : N+ — >N+ , such that for each integer 
n G , for all x G V^{= {z\\ z e V,%Ij{z) ^ 0}), alls g]0,£:o]- 



Ue[x) < (maxMe 



7(n)(C£:) 3 



It follows that at any xG V^, u^fx) tends to zero faster than any power of e and the 
convergence is uniform on any compact in V^. 

Proof. 

The proof will use be an induction on n. Given any smooth functions f, h on V we have: 



f[eAgh + ee{h)h] = sAgifh) + eemfh) - fh 

f[sAgh + ee{b)h] = eAgifh) + e9(b){fh) - fhe 

on V/, 



Ag/ , e(b)f 
f ~^ f 



(72) 
(73) 



where 6 = 6 + 2^, V/ = xG V, f{x) ^ 0}. For each nG Z+, set v„ = ilf^u^. Write: 



n+l 



t/-" {sAgU, + eB{h)u, + (t/- + v^c)m,] = ^—^e \AgU, + e{h)u,\ + (V- + v^c)m,V^", (74) 

on V^. 

Applying the formula ()72|) to ()74|) taking f='?/;"+^ and h=U£, we get, for n> 2: 
1 



AgV^+X + + Vn+\ + ^C^Vn - (n + \)e\d{^)i}) + Agl\)\Vn-\ 



-(n + l)(n + 2)£||V^||X~2 = 0. 



(75) 



Let M„ = max v„. Note that Mo=max'U£. Choose a point PgV such that M„+i=v„_(_i(P). 



Note that ■?/'(P)7^0. Evaluate the equation (f75|) at P. 

e \AgVn-,i{P) + ^(6)t;„+il (P) = £Ag^;„+i(P) > 0. 
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For n> 2: 



-{n + l){n + 2)e\\Vi^\\l{P)vn-2{P) < (76) 
The relation()76|l implies the inequahty: 



M„+i < VCeMn + {Ce){n + l)M„_i + {Ce){n + l)(n + 2)M„_2,for n > 2, (77) 
where : 



C = max( sup \ce\,snp[^^/\e{b)■^|J\ + \Agtp\,\\Vtp\\gi 

Vx[0,l] V 



The relation (fTTjl implies that: 



< V C7e- — + Ce ; h Ce- — , for n > 2 (78) 



(n + 2)! - (n + 1)! n! (n - 1)! 

aax(l, max?/', max-?/'^)Mo, 
remark and the reccurence relation (f7H|). 



Clearly M„ <max(l, max?/), max?/'^)Mo, if n=0, 1, 2. The lemma follows easily from this 
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